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Abstract 

This thesis comprises two parts centered around the functional renormalization- 
group framework: in the first part, I study the role of symmetries and conservation 
laws in approximate solutions, while in the second part I analyze Friedel oscillations 
and transport in Luttinger liquids with impurities. 

The functional renormalization group (fRG) has been developed as a new compu- 
tational tool in the theory of interacting Fermi systems. The effective behavior of a 
given microscopic model is calculated by solving coupled differential flow equations 
for the Green functions with an energy scale as the flow parameter. The symme- 
tries of the microscopic model imply Ward identities between Green and response 
functions. It is shown that solutions of truncated flow-equation hierarchies satisfy 
Ward identities if the cutoff bare action is gauge invariant. However, truncations are 
generally not self-consistent approximations in the sense of Baym and Kadanoff. 

The fRG is then applied to study Luttinger liquids. By computing the full spatial 
effective potential of a single impurity, long-range Friedel oscillations are observed in 
the density profile with the expected power laws for systems with up to 10 7 lattice 
sites. For a double barrier enclosing a dot region we find temperature regimes in 
which the conductance follows power laws with universal exponents, as well as non- 
universal crossover regimes in intermediate parameter regions. 

1 Version without the German summary; full version available online at 
jhttp : //elib . uni~ stuttga rt . de/opus/v olltexte/2005/2258/ 
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1 Introduction 



In one-dimensional metals electrons can move freely in one direction but are con- 
fined in the two transverse directions. The interaction between the electrons leads to 
Luttinger-liquid behavior with unusual properties different from conventional (Fermi- 
liquid) metals. In particular, the low-energy behavior of Luttinger liquids is strongly 
affected by impurities. Already a single static impurity has a dramatic effect: for 
a repulsive interaction, the backscattering amplitude grows as the energy scale is 
lowered, until at T = transport is inhibited and the chain is effectively cut into 
two pieces. The local density of states near an impurity, as well as the spatial density 
profile away from the impurity, obey characteristic power laws depending only on the 
bulk parameters. The conductance through a single impurity with varying parame- 
ters can be collapsed onto a single curve by a one-parameter scaling ansatz. A double 
barrier shows particularly rich behavior: it can be tuned to resonance, and additional 
scales are introduced by the separation of the two barriers and the detuning from 
the resonance. The conductance as a function of temperature is non-monotonous, 
exhibiting several distinct power laws, as well as a complex non-universal crossover 
behavior for intermediate parameter ranges. 

In recent years experiments on carbon nanotubes have allowed to measure the ef- 
fect of one or two impurities in an otherwise perfectly clean one-dimensional metal. 
While several field-theoretical predictions were confirmed, transport through a dou- 
ble barrier did not obey the expected asymptotic power laws. This led to a renewed 
theoretical interest to understand the behavior in intermediate parameter ranges 
accessible in experiments. Different analytical and computational methods applied 
to a spinless double-barrier model either supported or disagreed with the exper- 
imental data. This prompted us to investigate the problem with the functional 
renormalization-group method, which we have already used to treat complex multi- 
scale problems, such as Luttinger liquids with a single impurity. 

The functional renormalization group (fRG) has been developed in recent years 
as a new computational tool to study interacting Fermi systems. It is particularly 
efficient in low dimensions. Starting from a specific microscopic model, high-energy 
modes are successively integrated out to obtain the effective behavior on all energy 
scales. The method captures universal scaling laws in certain limits, as well as non- 
universal crossover phenomena at intermediate scales. 

Formally, the fRG flow equations constitute an infinite hierarchy of coupled dif- 
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1 Introduction 



ferential equations which describe the change of all Green functions as the energy 
scale is lowered. This hierarchy of flow equations produces the exact solution to all 
orders of perturbation theory at the end of the flow. In practice, however, the full 
hierarchy has to be truncated by neglecting the flow of some higher Green functions, 
which is justified perturbatively for weak renormalized interactions. In contrast to 
other renormalization-group methods we not only follow the flow of a few running 
couplings but of whole functions, such as the impurity potential. 

For the Id problems with one or two impurities, we approximate the interac- 
tion by an effective nearest-neighbor coupling but retain the full effective impurity 
potential. Our method is thus perturbative in the renormalized interaction but non- 
perturbative in the impurity strength. This already yields the expected universal 
scaling of the local density of states. In order to obtain the spatial density profile it 
is necessary to follow the flow of the density-response vertex. We treat the simpler 
case of spinless fermions; the more realistic modeling of electrons by spinful fermions 
is currently being considered. Our fRG results have been checked against numeri- 
cally exact density-matrix renormalization group (dmrg) data for systems with up 
to 1000 sites. The computation of the conductance as a function of temperature 
requires several extensions of the method. We develop the flow of the full impurity 
potential at fixed, finite temperature, as well as the temperature-flow scheme with 
self-energy feedback. Moreover, it is shown that on the level of our approximation, 
no corrections to the current vertex appear in the Kubo formula for the conductance. 
We are thus able to compute the conductance consistently within one approximation 
over several orders of magnitude in temperature, for arbitrary impurity strength. 

Besides the formal developments, the practical feasibility of the method depends 
crucially on the required computation time. Using a little-known mathematical 
theorem, an algorithm has been developed which scales linearly in the system size, 
instead of quadratically. For a lattice of 60.000 sites, the zero-temperature flow now 
takes minutes instead of days, and systems of up to 10 7 sites have been computed. 
This allows to find interesting regions in a large parameter space much more quickly. 

The vanishing of current-vertex corrections to the conductance is an example of 
a more general topic: the role of symmetries and conservation laws in the fRG for- 
malism. The microscopic model considered above, for instance, has a local U(l) 
gauge symmetry which implies charge conservation. As a consequence, the exact 
Green functions are related by Ward identities. In particular for transport calcula- 
tions it is crucial to respect these identities exactly even in approximate calculations. 
This raises the question whether typical approximations in the fRG, especially the 
truncation of the infinite hierarchy of flow equations, satisfy — or can be made to 
satisfy — Ward identities. This problem has been addressed ten years ago in the 
context of gauge theories with a fluctuating gauge field, using either modified Ward 
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identities or the background-field method, but simple gauge-invariant truncations re- 
mained elusive. In those cases where a gauge-invariant flow is possible, for instance 
if the model is regularized not by a momentum cutoff but by a finite temperature, 
we show that even truncated flows satisfy the Ward identities. On the other hand, 
we find that "self-consistency" between Green functions of different degree, a feature 
of the conserving approximations by Baym and Kadanoff, is generally not satisfied 
by common truncations of the fRG flow equations. 

This thesis is organized as follows: 

• In Chapter 121 the fRG formalism is introduced. After a brief review of generat- 
ing functionals and their expansion in terms of Green or vertex functions, an 
infrared cutoff is defined which introduces a scale dependence in the generat- 
ing functionals. A derivative with respect to this scale leads to functional flow 
equations. These are then expanded in terms of their constituent Green func- 
tions to obtain an infinite hierarchy of coupled differential flow equations for 
the Green functions and a diagrammatic representation of the flow equations. 
The merits of different schemes are compared. 

• In Chapter El Ward identities are derived expressing the symmetry of the bare 
action in the functional formalism. A momentum cutoff generally modifies the 
Ward identities. For other flow schemes which preserve Ward identities we 
show that they hold even in truncated flows. Conserving approximations are 
reviewed as an example of self-consistent approximations. It is then shown 
that common truncated fRG flow equations are generally not self-consistent. 

• In Chapter 0] the general fRG formalism is applied to study one-dimensional 
correlated fermion systems (Luttinger liquids) with impurities, in particular 
their single-particle and transport properties. In this technical part the precise 
form of the flow equations on the lattice is derived, as well as the details 
of the finite-temperature cutoff procedure, truncations of the flow-equation 
hierarchy and parametrizations of the flowing vertices. At the end of the flow 
we obtain the effective impurity potential (self energy) and the renormalized 
density profile. In order to compute transport in our approximation, we then 
have to solve the scattering problem of non-interacting electrons in this effective 
potential. It is also shown that current-vertex corrections to the conductance 
vanish in our approximation, in accordance with the Ward identities. The loop 
algorithm for nearest-neighbor interaction, which scales linearly in the system 
size, is derived in appendix iBl 

• In Chapter El new results are reported for the Friedel oscillations of the spatial 
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density profile generated by a boundary or impurity in one dimension, and the 
temperature-dependent conductance through a double barrier. For appropri- 
ate parameter ranges universal scaling is observed with several distinct power 
laws in temperature. In intermediate regions the full non-universal crossover 
behavior is obtained, suggesting an interpretation of recent measurements on 
carbon nanotubes. 

The publications based on this thesis are listed on page 11171 
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2 Functional renormalization group 



Challenging many-body problems often involve effects on many energy scales. In 
perturbation theory one has to perform loop integrals over all energy scales, which 
may lead to infrared or ultraviolet divergences. Some of these divergences have a 
physical origin indicating for instance a phase transition, while others are an artefact 
of perturbation theory. Wilson's exact renormalization group (rg) |Wilson 1971[ 
WilsonfcKogut 1974| provides a method to deal with such problems: the different 



energy scales are successively taken into account by integrating out momentum shells. 
This can be done by introducing for example an infrared cutoff in the bare propagator 
which suppresses all modes with an energy below the cutoff scale A. Then, all 
correlation functions depend on the scale A. One follows the change (flow) of the 
correlation functions as the cutoff scale is lowered until finally the cutoff is removed 
and the original theory is recovered. An important advantage of this procedure is that 
the right-hand side (rhs) of the flow equation remains regular even if perturbation 
theory leads to unphysical divergences. 

There are several variants of the exact functional RG flow equations. After Wil- 
son's early review WilsonfcKogut 1974| on the exact RG, |Polchinski 1984] derived 
equivalent continuum flow equations with a smooth cutoff in order to prove per- 
turbative renormalizability of massive Euclidean tp 4 theory in D = 4. Keller, Kop- 
per, and Salmhofer |kks 199*2] simplified and extended the proof and showed that 
these flow equations determine the connected amputated Green functions (cf. sec- 
tion I2.2.4j) . |WegnerfcHoughton 1973] derived flow equations for a sharp cutoff, 
however there were ambiguities which they avoided by assuming discrete momenta. 
|Weinberg 1976 took the continuum limit and found that the flow could be formu- 
lated without ambiguity by expanding the connected amputated Green functions 
in trees, i.e., in terms of one-particle irreducible (Ipi) vertex functions. The Ipi 
flow equation of the Legendre effective action, the generating functional of the Ipi 
vertex functions, was derived by |NicollfcChang 1977} IWetterich 1993] Ibdm 1993] 
IMorris 1994] SalmhoferfcHonerkamp 2001] KopietzfcBusche 200 1| . Another scheme, 
obtained from the Polchinski scheme by Wick-ordering jWieczerkowski 1988]ISalmhofer 1998] 
ISalmhofer 1999] . is particularly suited for rigorous proofs because it allows strong 
bounds on the growth of correlation functions, even near the Fermi surface and to 
all orders in the renormalized interaction. 

By then, several people had started to use the functional RG (fRG) for fermionic 
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lattice models. This introduces several problems not present in ip 4 theory, for exam- 
ple the determination of the Fermi surface, whose shape and position is not known a 
priori. Another problem is that even at low energies, when only momenta close to the 
Fermi surface are important, the two-electron interaction is a complicated function of 
momenta. Therefore, it needs to be parametrized by many discrete couplings, in con- 
trast to the single renormalized coupling A at zero external momenta for the (p A the- 
ory. On the other hand, the lattice provides a natural ultraviolet cutoff which leads to 
significant simplifications as compared to the continuum field-theoretical models. Im- 
portant applications of the fRG in condensed-matter physics include the 2d Hubbard 
model using the Polchinski scheme jZanchifcSchulz 1998[ IZanchifcSchulz 2000| . the 
Wick-ordered scheme |HalbothfcMetzner 2000j and also the Ipi scheme |hsfr 2U0*T] . 
while Id impurity problems and Luttinger-liquid physics are conveniently investi- 
gated in the Ipi scheme jMMSS 200"2~al Immss 2002bt Iaemmss 2U04j . 

Initially, many of the flow-equation schemes have not been derived in the most 
straightforward way, and I found it worthwhile to derive them again in a simple 
and uniform notation, highlighting the relation between different schemes. While 
the results are not new, some derivations are much easier than those found in the 
literature, and I hope the reader new to this method will find them helpful. After 
introducing the functional formalism and several types of correlation functions in 
section 12.14 the most commonly used flow equations are derived in section 12.21 I 
compare important features of the different schemes in section |2~31 

2.1 Functional formalism 

2.1.1 Bare action 

A system of interacting spinless 1 fermions is described by the action 



A 

where for a translationally invariant system (|2.3|) . the multi-index K = (u,k) con- 
tains space-time indices (frequency, momentum) and could be extended by internal 

1 Ultimately we aim to describe electrons which are fermions with spin; presently in the applications 
in Chapter^] however, we only consider fermions without spin and, therefore, specialize to this 
case. 



with the kinetic (quadratic) term defined as a scalar product 






(2.3) 
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2. 1 Functional formalism 



degrees of freedom like the spin projection a. The ipK, 4>k are Grassmann variables, 
and Q(K) = iuj — £ fc is the inverse bare propagator with £ fc = e k — \x the dispersion 
around the chemical potential \x. If, however, one includes an impurity potential in 
the bare propagator that depends specifically on space and not just differences of po- 
sitions, translational invariance is broken: then one has to use the general expression 
(|2.2|) . cf. section l4~Tl The bare propagator is the inverse operator of Q, 

C:=Q-\ (2.4) 

The functional Vo[ip, ip] is the bare many-body interaction, for example in the typical 
case of a density-density interaction (with frequency conservation implicit): 

V [i>,ij] = / dxdyV (x - y) n(x)n(y) (2.5) 



kx,k 2 ,q 

2.1.2 Partition function and connected Green functions 

All information about the physical system with action (|2.1|) is encoded in the nor- 
malized partition function 

Z[r), f)\:=^r [ [ditf>] ^ M e^H^ , (2.6) 
Zq J 

a functional integral with weight e s and coupled to Grassmann source fields r)K, 
f)K- The integration measure is abbreviated as [diptp] := Y\ K dipx dtpx- The parti- 
tion function Z[rj, fj] is the generating functional for the Green functions (connected 
and disconnected). The normalization constant Z is the non- interacting partition 
function, 

Z ■= I [d^] ^k^kQk^k = y f dij K d^ K e^ KQK ^ K = J[Q K = det(Q). 

J K J K 

(2.7) 

It is convenient to absorb the quadratic part of the action as well as the normalization 
factor into the measure, 

Z[ v , fj]= J dn Q [i/;, $ e~ Vo ^ e^H^) (2.8) 

with the normalized Gaussian path-integral measure 

dfi Q [i/j,ij>] := i#$e(^ w ) such that / dfi Q [i/j,^} = 1. (2.9) 

Zq J 
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If there is no interaction, Vq = 0, the integral is quadratic (Gaussian) and can be 
performed analytically by completing the square: 

Z nonint hv} = J dfi Q [^] e -6M-(^) 

= — / [diH>] e (*W)-(^)-(^) 

= — [ [dtptp] e (li>- ctf ilQl4>-Cv})-(fi,Cv) 
Zq J 

= e-^ c n) J diigty -Cri,$- 0% 
= e -(v,cv) 

where C* is the transposed propagator. 

The connected Green functions are generated by the functional 

£[77,77] := - In £[77,77] 

as 

00 ^ m 

£fo>*j] = J^T^fia Yl Gm{K' 1 ,...,K' m ,K 1 ,...,K m )Y[f] K/3 r] Kj 

m=0 ^ V K!...K m K[...K' m j=l 

and 

G m (K[, . . . , K' m , Kx,..., K m ) = -(lp K [ ■ ■ ■ ^K' m ^K m ■ ■ ■ did) cam 

dr] Kl . . . dr] Km 6r) K!m . . . br)K[ v=v=o 
respectively. For the non-interacting system follows 

0**** = (fj, Crj) (2.10) 
such that Gf^^^K) = C(K) is the bare propagator, and all other G£° nmt vanish. 



2.1.3 1PI vertex functions 

There is another set of correlation functions that is particularly useful for describing 
phase transitions and fields whose expectation value does not always vanish, such as 
the order parameter in a symmetry-broken phase: the one-particle irreducible (Ipi) 
vertex functions •y m generated by the functional T[cj), </>]. T is obtained via Legendre 
transformation from the connected Green functions |Zinn- Justin 2002{ section 7.8], 

T[0,0] + (0,Q0) := £[77,77] + (0,7/) - (77,0) . (2.11) 
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This differs from the textbook definition by taking the inverse bare propagator Q 
out of the one-particle vertex function. The transformations between and r\ are 



SfT - Q*<l> = fj, 5 v g = (f>, S^G = (SfyT + Q)- 1 
S^T + Q(f) = ri, 5$ = (f). 

T is expanded in <j) K , (fi K as 



(2.12) 



oo - m 

r [^0] = 5Z?^nI Yl Yl lrn,{K[,...,K' m ;K 1 ,...,K m )\^4> K ^ Kj . 

m=0 \ K 1 ...K m K[...K' m j=l 

Each Ipi vertex function 7 m is made up of those diagrams of G m which cannot 
be split into two disconnected parts by cutting a single line. As all self-energy 
contributions on external legs are one-particle reducible with respect to the main part 
of the diagram, full propagators are amputated from all external legs. In the special 
case without interaction with Q nonmt [r],f]] = (f},Crj) we obtain r nonmt [0, </>] = 0. 

2.1.4 Connected amputated Green functions 

Another way of looking at a system is by considering the generating functional of 
the connected amputated Green functions, the effective interaction V[x, x]'- 



e -v[ x ,x] := J d/jiQ[M] e -W+>^+x] (2.13) 



<fi=ip=Q 



e -Vo[8p,6 v ] e &,Cip) e (<p,x)-(x,V>) 



e ~V [50,6^,] e (S x ,C6x) e (v,x)-(x,v) 



tp=tp=0 
ip=if=0 



e Ac e ~ Vo[x '^ 



where the functional Laplace operator is defined as 

A ( ^ C ^ \ f d d ^ C( ) ^ ^ C ^ 

°' \$X &XJ J Sx{y) ' Sx{x) k ^ Xk K 

I shall use the shorthand notation only if it is unambiguous on which Grassmann 
variable the derivatives act and use the explicit notation otherwise. Ac acts on a 
functional J-[x, x] m the following way: the derivatives j- and ^ pick an ingoing 
and outgoing leg from each diagram in T and connect them by a bare propagator 
C. In this picture, perturbation theory for e~ v may be visualized as follows: e~ Vo 
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is a collection of any number of disconnected bare interaction vertices, and Ac e~ Vo 
contains all diagrams where either Vq is closed by one C loop to create a tadpole 
diagram, or two Vo's are connected by a C line to form a tree diagram. Repeating 
this procedure indefinitely, higher tree and loop diagrams appear, and finally e~ v = 
e Ac e~ Vo contains all Feynman diagrams. Taking the logarithm to obtain V, only 
the connected diagrams are retained |Zinn- Justin 2002[ section 1.2.1]. 

What is the relation between Q and V as both generate all connected Green func- 
tions? We observe that by the substitution x '■= Crj, % '■= C l fj, 



e -v[c v ,c% = I d ^ ) $] e -vM+Cv,*+c'fi] 

duqW - Crj, $ - C% e ~v*Wm 



— e (v,Crj)-Glv,v] 

such that 

V[Cr ] ,C t r ] } = g[r ] ,r ] }-(r ] ,Cr ] ). (2.14) 

In the non-interacting case, V = because (fj, Crj) cancels the non-interacting part 
(I2.1()j) of Q, in accordance with the name effective interaction. Generally, V generates 
connected amputated Green functions V m , 

oo m 

V ^^ = Hj-^ E E V m (K' 1 ,...,K' m ,K 1 ,...,K m )HxK' j XK j , 

m=0^ K 1 ...K m K[...K' m j=l 

(2.15) 

i.e., connected Green functions G m with bare propagators C(K) amputated from all 
external legs: 



V m (K[,...,K' m] K u ...,K m ) = - — - V[x,x] 

OXKi ■ --0XK m OXK' m ■ ■ -0Xk{ 



x=x=o 



GmiK'i, . . . , K' m ) Ki, . . . , K m ) S m: i 



C t {K[)---C{K 1 ) C(i^i)' 
2.2 Renormalization-group flow equations 

When computing a Green function perturbatively by summing the contributions of 
certain Feynman diagrams, there can be infrared divergences both as an artefact of 



11 



2.2 Renormalization- group flow equations 



perturbation theory and physically as an indication of a phase transition. These can 
be regularized by an infrared cutoff A in the bare propagator which suppresses modes 
of low frequency or momentum close to the Fermi surface. The change of the Green 
functions as the cutoff scale is changed is governed by the renormalization- group flow 
equation. 

Alternatively, one can regularize the problem by going to sufficiently high temper- 
atures T > T c , weak coupling g 2 U < U, and/or a finite system size N; then one can 
consider the flow of the Green functions as the parameters T or g are changed. 



2.2.1 Regularization and flow parameters 

The flow schemes used for the functional RG are constructed by making only the 
quadratic part of the bare action depend on the cutoff scale or flow parameter. In 
all cases the regularization is done with respect to energy scales. 



Frequency and momentum cutoff 

The frequency cutoff is defined by multiplying the bare propagator in the action 
with a cutoff function x A ( UJ ), 

C\K) = x A (uj)C(K), (2.16) 

where x A (o;) cu ^s out modes with frequency \cu\ < A in the frequency basis where 
C(K) is diagonal. It may be either a sharp cutoff function using the step function 

x A M = e(M-A), 

or a smooth cutoff. Generally one would think that a smooth cutoff with a differen- 
tiable cutoff function is easier to handle. However, especially at zero temperature, a 
sharp cutoff can in fact be favorable (cf. section 12. 2. 3|) . 
The momentum cutoff is defined analogously with 

X A (fc) = 6(16*1 -A). (2.17) 

Generally, the propagator is split as 

C(K) = C A {K) + D A {K) (2.18) 

where C A (K) is the propagator for the high-energy (hard) modes, and D A (K) is the 
propagator for the remaining low-energy (soft) modes that are yet to be integrated 
out. As A — > 0, the cutoff is removed and C A — ■> C, D A — > 0. 
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It is simplest to define the cutoff in the basis where C is diagonal, as a function 
multiplying the kinetic energy eigenvalues e with x A ( e )- We will later (section l3.2.2|) 
see an example where in the presence of an external field the propagator is not 
diagonal in momentum space: then the cutoff can be defined by diagonalizing first. 

Note that at finite A, the frequency cutoff leads to non-analyticities in certain 
correlation functions, for example in the current- current response function. This 
disadvantage is absent in the temperature flow. 

Temperature flow 

In order to use the temperature T as a flow parameter in the fRG formalism, only 
the quadratic part of the action may depend on T. Initially, however, the action has 
powers of T both in the kinetic and in the interaction part: 

iu n k 

+ 2 T3 E V^K^K^k^k 

K lt K 2 ,K[ 

By rescaling the fields tp := T 3 / 4 ?/>, one can shift the explicit T dependence of a 
quartic interaction into the quadratic part |HonerkampfcSalmhofer 2001| : 

tUJn k 

+ 1; Y V ( K 1> Kl > K *) '4'K' 1 1pK 1 +K 2 -K' 1 1pK 2 i'K 1 ■ 
Kx,K 2 ,K' x 

We can now perform the fRG flow on this action to obtain Green functions in 
terms of the rescaled fields ip for a whole temperature range and rescale them back 
to obtain the usual Green functions 2 

G m {...)=T^I 2 G T m [...). (2.19) 
Interaction flow 

There are systems which can be treated also by a much simpler method |hrae 2UU4j . 
namely by simply rescaling the propagator with a number g G [0, 1]. Then, all 
Feynman diagrams become dressed by a global factor of g for each internal line. This 
seemingly trivial change to weight each diagram by a power of g has the advantage 

2 This is in contrast to |HonerkampfeSalmhofer 2001| where the powers of T in the original action 
are considered part of the Green functions. 
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that the fRG flow equation for the Green functions from g = to g = 1 resums 
infinite subclasses of Feynman diagrams which would otherwise be more tedious to 
do. At the same time we can rescale the fields and observe that for g < 1 the model is 
the same as one with full propagators but reduced interaction strength g 2 \U\ < \U\: 

S~($,[Q/g\ij,)-UWiH, 
~ y^,Q'ip \ — g 2 U Tp^ipijj 

where ip = g -1 ! 2 ^. Integrating from g = to g = 1 thus yields a flow in the bare 
interaction. 



2.2.2 Connected Green function flow 

The above regularization and flow schemes lead to an action with a quadratic part 
depending on a scale parameter A, which for notational simplicity shall include T 
for the temperature flow and g for the interaction flow, 

S A [^] := ('0, Q A ip) - V [i>,i>] , (2.20) 

where the inverse bare propagator with cutoff is defined in analogy to equation 1)2.40 . 

Q A : = (C A )-\ (2.21) 

Generally, if x A (K) = 1 then C A (K) = x A (K)C(K) is just the original propagator 
while for x A = 0, C (K) = and the kinetic term (C A ) _1, 0) in the action 
becomes infinite, giving the cutoff modes of the fermion fields an infinite mass such 
that they are frozen out. By the action S A all generating functionals, and thus their 
Green function components, depend on A. The generating functional Q K for the 
connected Green functions in the presence of the cutoff is defined as 

e 



-G A lv,H}. = I ^ QA [^^] e -^[^] e -(^)-(^) (2.22) 

0,qH)-v o [iPM e -($,v)-{v,il>) (2.23) 



7 A 
z 



where Q is replaced by Q A , and the normalization factor is changed accordingly to 
Zq = det(Q A )- The flow equation for e~ gA is obtained by taking the A derivative, 
denoted by the dot, on both sides of equation ()2.23jh 



(d A G A ) e- gA ^ = - ^7^ A e- gA + 

Lit- L \cJ 



+ 



d\ det Q _ g t 
letQ A 6 
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(-Tr(g A C A )- A QA )e^ 



where the first term comes from the derivative of the normalization factor, d\ In det(Q A ) 
Tr dji ln(g A ) = Tr(g A C A ). Here, Tr denotes a sum over all space-time indices. Thus, 
the flow of Q A is 



d A G A [ V , fj] = Tr(Q A C A ) — Tr g , _ 11 + If-i^l, Q 



SrjSr] J \ §7] 5rj 



(2.24) 

As a check, in the non-interacting case 

g A [r ] ,r ] } = (r ] ,C A r ] ) 
r) 2 G A 

q\ _ y _ q\ 

1 5r] 5fj 

G A = Tr(Q A C A ) - Tr(Q A C A ) + (-fj, C A Q A C A V ) = (fj, 6%) . 

Two technical notes are in order. Naively, the Q A appearing in the flow equation 
flXHD for Q A looks ill-defined: 

to - -Smp QiK) 

contains a division by zero for the cutoff modes where x A (i^) = 0. But all Green 
functions G A have C A on their external legs so only the combination 

S A {K) := C A Q A C A = -C A = -x A {K)C{K) {single-scale propagator) 

(2.25) 

appears. This quantity is well-defined and, for a sharp cutoff, has support only on 
the scale A since x A is a 5 function. The other seemingly ill-defined contribution is 

Tr(g A C A ) = -X:^||, (2.26) 
but this is canceled by the bare part of the second term in the flow equation ([2.24)1 . 

5 2 g A _ 

77=77=0 



Tr(Q A C A ) - Tr (q a = Tt(Q a [C a - G A }) 

V Or] Or]/ 77=77=0 

= -^Tr(g A C A [S A C A ] fe ) 



k=l 

which again contains only the well-defined combination S A = C A Q A C A . On the 
second line, the full propagator G A has been expanded into a geometric series via 
the Dyson equation (G A )^ = (C A )~ 1 - S A . 
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2.2.3 1PI vertex function flow 



From the flow equation (|2.24jl of the connected Green functions it is simple to derive 
the flow of the Ipi generating functional T A . We use the Legendre transformation 
(|2.12j) but all functionals and Q are taken at scale A, 



r A [0,0] + (0,Q A 0) := g A [r),fj] + (0,77) - [f)A] 
We express Q A and its derivatives by T A , 

<9 A r A [0,0] = d A g A [r),fj}- (0,qV) 

,Q^-^)-TrQ A 



(2.27) 



W A hv] AA SG A hv}\ ^ AA [S*g A [ri,fj} 



5r) 



dij 



5 2 r A \ 



Q 



5rj 5rj 
-C A 



C 1 



Q 



•SA, 



Q 



to arrive at the lPI flow equation 



d A T A = TtQ a 



5 2 T A 



5<p5<p 



Q 1 



(2.28) 



If we define the abbreviation T^> := S 2 T A 
in a more compact form as 



the flow equation can be written 



<9 A r A = d A Tr In Q A - d A Tr ln(Q A + T^) = d A Tr In qJ^® 



(2.29) 



where in equation ()2.29j) . the A derivative acts only on the cutoff functions inside 

Q A - 

The inverse of the functional {Q A + T^) is the inverse full propagator in the 
presence of the fields 0, and can be safely defined via a geometric series, as in 
the Dyson equation. To see this, split (Q A + r^ 2 ^) into a part independent of 0, 
which is the usual inverse full cutoff propagator [G A ] _1 = Q A — S A , and a remaining 
functional r A f0, 01, 



2rAr 



Q 1 



6 2 T 



A\-l 



+ r 



yA _ A_a 2 r A [0,0] 

-Zu — 7i — " 



0=0=0 



Thus, 



2rAr 



Q 1 



6 2 T 



5050 



(1 + G A f A ) _1 G A = ( 1 - G A f A + [G A f A ] 2 )G / 
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Together with the definition of the single-scale propagator S A in analogy to (|2.25jl . 

S A := G A Q A G A (single-scale propagator), 
the flow equation becomes 

<9 A r A = Tr(Q A [C A -G A }) + Tr(S A [f A — f A G A f A +■■■]) (2.30) 

where the last trace consists of a one-loop term with any number of vertices T A , which 
contribute at least two external legs each, connected by one single-scale propagator 
S A and several full propagators G A . Writing the flow equation in terms of the 
components 7 A , we obtain the following diagrams for the general hierarchy of flow 
equations and for the first few levels of this infinite hierarchy: 
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In the general form of equation (|2.31[) . the wiggly line in the last diagram denotes that 
the expansion of the inverse second derivative yields further 1-loop diagrams with 
an appropriate number of higher vertices insertions, such that the number of exter- 
nal legs on the lhs and RHS matches. This flow scheme is derived in |Weinberg 1976[ 
NicollfeChang 19771 IWetterich 19931IBDM 19931 [Morris 1994(|SalmholerfeHonerkamp 2001| . 



Sharp-cutoff flow equations 

Consider the flow equation (J2.29J) in a basis where Q is diagonal and the cutoff 
function x A is multiplicative so that it can be taken out of Q A , 

Assume that all x A along the loop have the same value. This is clearly not the general 
case but is sufficient in our application in Chapter 0] where x A is a sharp frequency 
cutoff and all vertices are evaluated at zero frequencies, such that all propagators in 
the loop have the same frequency and hence the same value of the cutoff function. 
Then <9aX A = — <KM ~~ A) restricts the frequency to the A shell but the RHS contains 
0(|o;| — A) which has a step right at the shell. Such expressions are unique if the sharp 
cutoff is implemented as the limit of increasingly sharp, broadened cutoff functions 
e with broadening parameter e — > 0. This is demonstrated by a lemma due to 
Morris |Morris 1994] : for an arbitrary continuous function f(t), 

5 t {x - A) f(Q e (x - A)) -> 5 e (x - A) T f(t) dt (2.33) 

Jo 

where 5 t (x) = Q' e (x). Then, equation (|2.32|) reads 
<9 A r A = Trx A C dt^hi Q 



dt q + tr< 2 ) 



Tr x A l n 



Q + r( 2 ) 



rp -A, Q Q~Z A 

lr x In 



Q - s A q - s A + r A 

Tr - X A [ln(l - C£ A ) + ln(l + G A f A 
i-^tr[ln(l-C S A )-f:^((5 A fY 

ui=±A k=l 



(2.34) 
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where the T = Matsubara sum ^ J dui has been performed in the last line such 
that tr denotes the sum over the remaining spatial indices. The new sharp-cutoff 
propagator is defined as 

G A := [Q - E A }-\ (2.35) 

This propagator has no step at = A as opposed to G A , hence the RHS of the 
flow equation which contains one-loop terms built from powers of G A r A is a smooth 
function of A and the vertices T A . 

However, the sharp frequency cutoff is only possible if the loop integral and hence 
the vertices T A are continuous: at T > 0, an imaginary frequency integral is restricted 
to discrete Matsubara frequencies uo = uo n at which all vertices have a step as a 
function of A, hence the condition that f(t) be a continuous function in Morris' 
lemma ()2.33j) is not satisfied any more. 



2.2.4 Connected amputated Green function flow 

The Polchinski scheme |Polchinski 1984j has been introduced to prove the renormal- 
izability of the (p A theory via the flow of the effective interaction, the generating 
functional of the connected amputated Green functions (cf. section l2.1.4j) . In the 
presence of a cutoff, the effective interaction V A [x, x] is defined as 

e -v A [ x ,x] .-J dfiQA^] e -W+x,^+x] 

= e (£' cA *) e~ v ^ (2 ' 36) 
= e A c A e~ Vo ^ x '^ . 

Thus, all Feynman diagrams contain not C but C A on the internal lines. In order to 
obtain the full V without cutoff, one has to apply the missing e A o A on e~ yA ; 

e -V[x,x] = gAc e -Vo[ X ,x] 

— g A DA + C A e ~ V olX,X] 

= e A D A e A c A e - Vobc >% (2.37) 

= e A DA e -V A [ X ,x] _ 

Equation (J2.37)) can be interpreted as follows: if all Feynman diagrams (connected 
and disconnected) with hard internal lines are reconnected again in all possible ways 
with soft lines, one obtains all diagrams with full lines. 

The flow equation for V A is derived by taking the A derivative of equation ()2.36|) . 

a -\ )A V A a -V A 
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= _ e V A 9A(e A cA e -V 0) 

= -e vA A cA e~ v \ 
to arrive at the Polchinski flow equation 



«9 A V A = Trfc A i-4) 



$x ' &x 



(2.38) 



In a more compact notation, 



d A V A = A cA V A - \ A^ A V A [ X i, Xi] V A [X2, H 



(2.39) 



Xl=X2=X,Xl=X2=X 



where A^ 2 = C~s^j + (dfj'^fr) connec ts two different vertices. The initial 
condition is the bare interaction: 

V Ao [x, x] — ^o[Xj x] — y (f° r a two-particle interaction). 

The graphical representation of the Polchinski equation for the connected amputated 
Green functions V A features both tadpole and tree diagrams: 






C' 




V A 

v m-k+l 



2.2.5 Wick-ordered Green function flow 

The connected amputated Green functions are the expansion coefficients of the gen- 
erating functional V in terms of monomials of the source fields x-, X-, whereas the 
Wick- ordered Green functions are the expansion coefficients in terms of Wick-ordered 
polynomials of the source fields, e A o A (^^) m . The construction is completely anal- 
ogous to the use of Hermite polynomials e d *x m as compared to the monomials x m : 
any analytical function f(x) can be expanded uniquely in terms of x m (Taylor ex- 
pansion) or in terms of Hermite polynomials, which provide another complete or- 
thogonal basis. In appendix El I solve the heat equation in real space using Hermite 
polynomials; while this is different from the usual textbook solution, it provides 
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a low- dimensional and intuitive example of the functional formalism. The Wick- 
ordered generating functional is defined as |Wieczerkowski 19881 ISalmhofer 19981 
ISalmhofer 19991 SalmhoferfeHonerkamp 2001| 



W A [x,x] = e A - A V A [ X ,X]- 



(2.40) 



Diagrammatically, the effect of Wick ordering is that the coefficients W A of W A 
contain all possible D A loops (tadpoles) on the connected amputated Green functions 

yA. 






w 



A 
m 



V A 

v m+l 



V A 

v m+2 



The flow equation for W A can be derived from the Polchinski flow equation ()2.39J) 
with the replacement A^a = — A^ A (since C A + D A = C is independent of A), 

8 A W A = (d A e A ° A ) V A + e A « A d A V A 

= A DA e A ^V A + e A ^(-A DA V A ) + e A ^(iA^ 2 A V A V A ) 
= exp{A^ 12+22 }|Ag A V A V 2 A 
= i(Ag A e A ^)(e A ^V A ) (e A ^V 2 A ). 

The first two terms in the second line cancel each other, and A^ 12 " 1 " 22 = A]^ A + 
A^a+A^a contains derivatives with respect to the fields labeled 1 and 2, respectively. 
Finally, we obtain the Wick- ordered flow equation 



±<9 A (e A ^)W A W 2 A 



(2.41) 



The initial condition is 



W A( '[x,x] =e Ac V [x,x} = g 



^ (for a two-particle interaction). 



The Wick-ordered functional W A differs from V A only for A > 0: as A — > 0, W A — > 
W = V, i.e., W A ° has a different starting point (initial condition) from V A ° but flows 
via a different route to the same final functional (cf. figure l2~T]) . The diagrammatic 
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representation of the flow equation for the components of W is 



i) A 




Note that all terms on the RHS of the flow equation are bilinear in the vertices 
and are either tree or higher loop diagrams. In the exact hierarchy without trun- 
cation, the flow equation of each vertex has infinitely many terms on the RHS with 
higher vertices connected by many loops. One internal line is the single-scale prop- 
agator D A , while all others are soft-mode propagators D k . As the cutoff scale is 
lowered, only momenta in a small neighborhood of the Fermi surface appear on the 
internal lines. This justifies a particularly efficient parametrization of the coupling 
functions which depend on momenta anywhere in the Brillouin zone: at a low cut- 
off scale, all internal momenta are close to the Fermi surface, so it is sufficient to 
parametrize the couplings by their values with all momenta projected onto the Fermi 
surface. This parametrization is employed in computations of the 2d Hubbard model 
IHalbothfcMetzner 20001 iRoTiefcMetzTier 2005]. 

2.3 Summary 

In this chapter I have derived the functional renormalization-group flow equations 
in the Ipi, Polchinski and Wick-ordered schemes. While this is not new, I believe 
that the derivations are formally simpler and more straightforward than in much of 
the literature. Furthermore, the treatment of the sharp cutoff on the functional level 
(cf. section 12. 2. Hj) has not yet been published to my knowledge. 

Figure 12.11 provides an overview over the different fRG schemes: in the center 
is the Polchinski scheme of connected amputated Green functions V A which have 
essentially the same structure and flow equations as the connected Green functions 
Q A . By Legendre transformation to T A we have obtained the Ipi scheme which 
starts from the same initial condition but parametrizes the physical properties in 
a different way particularly suited for symmetry breaking. It has the advantage 
that the internal lines are full propagators taking into account all self-energy effects 
already known at scale A, while it has the disadvantage that one has to integrate 
the internal loop over the whole Brillouin zone. 

On the other side, the Polchinski scheme leads to the Wick- ordered scheme W A 
which starts from the first-order Hartree-Fock solution as the initial value of the self 
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Figure 2.1: Overview over the relation between the different flow- equation schemes. 

energy but converges towards the same connected amputated Green functions V A at 
the end of the flow. Its main advantages are that internal lines are restricted to a 
small neighborhood of the Fermi surface, and that its simple power counting allows 
rigorous proofs of renormalizability to all orders in perturbation theory. 

The next difference arises with truncations. While the full hierarchy of flow equa- 
tions in any scheme leads to the correct solution to all orders in perturbation theory, 
the lowest orders in different basis sets of correlation functions capture different as- 
pects of the the solution. As an illustration, imagine a complicated real analytical 
function f(x), expand it in different sets of orthogonal polynomials and retain only 
the first few coefficients: the approximate result will look different for each basis set. 
A third question is whether the basis functions allow for efficient parametrizations in 
terms of a small set of coefficients, for instance by projecting onto the Fermi surface. 

Hence there is not one single flow scheme superior for all applications, but the 
method of choice depends on the particularities of the model studied. 
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In constructing our microscopic model we demand that it satisfies certain physical 
properties, such as charge conservation. This particular property is guaranteed by 
the requirement of U(l) gauge invariance: by the Noether theorem every continuous 
symmetry has a conserved quantity associated with it, in this case the electrical 
charge. In a classical calculation, the conservation of charge is expressed by the 
continuity equation 

d t p{x) + V- j{x) =d^{x) = 

where p(x) and j(x) are the charge and current densities, respectively. In quantum 
theory the current j^{x) becomes an operator and the continuity equation 

is evaluated inside an expectation value together with insertions of other opera- 
tors. If the time derivative acts not only on the current operator but also on the 
time ordering in the expectation value, it generates differences of Green functions 
|ItzyksonfeZuber 1980| . In momentum space with momentum transfer q to the ex- 
ternal field, and in units such that the electron charge e = 1, 




This is an example of a Ward identity (wi): these are constraints on Green functions 
expressing the underlying symmetry |Ward 1950 . For the one-particle Green func- 



tions equation ()3.1|) can be easily derived by hand, but it will become more tedious 
as more external legs are added, each corresponding to an operator insertion in the 
expectation value. Instead, we shall employ the formalism introduced in Chapter El 
to derive a functional Ward identity applicable to all m-particle Green functions si- 
multaneously. This will allow us to investigate general properties of Ward identities 
to all orders in the fields and in the interaction. 

Let me mention one important point already now in order to avoid confusion 
(the issue will be explained in more detail below): there are two versions of the 
Ward identities in the literature, one relating the difference of Green functions to 
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the current as above (response-function Ward identity), the other relating it to the 
two-particle interaction with a loop closed (self- consistent Ward identity). Both are 
equivalent in the exact solution and in conserving approximations but generally differ 
in other approximations such as truncated fRG flows. I will show that in truncated 
flows the response-function Ward identities can still be satisfied in the absence of a 
momentum cutoff, for instance in the temperature-flow scheme, while self-consistency 
and thereby the self-consistent Ward identities are in general violated. 

We proceed as follows: after the functional derivation of Ward identities in section 
13.11 we shall see how a momentum cutoff breaks response-function Ward identities 
and how this can be fixed in flow schemes without momentum cutoffs (section l3.2|l . 
To answer the related question of self-consistency (section 13.3)1 we first review the 
construction of self-consistent conserving approximations of Baym and Kadanoff. We 
then show that the self-consistent Ward identities are generally violated by truncated 
flow equations. The results are summarized in section l3~H 



3.1 Gauge invariance and Ward identities 

We couple the fermionic action to an external field for two reasons: it allows us to 
compute linear response by derivatives with respect to the external field, and it is 
necessary to make the action satisfy the desired symmetry. 

In order to compute electrical transport properties, we need to study the response 
of the system to the external electromagnetic potential. This is done by including 
in the action a term coupling a current operator j^(x) to the external field A^x), 

S[il),$,A]:= S[^$\-(j»,A») (3.2) 

where the summation over [i is understood. Then, the response functions are deriva- 
tives of the Green functions with respect to A fl (x). The current operator j^ty, ip, A] 
is a composite operator of the fields ip, ip and generally depends also on the external 
field A. We shall now demand a symmetry and construct the current such that this 
symmetry holds. 



3.1.1 Gauge transformation 

The local U(l) gauge transformation is |Zinn- Justin 20021 chapter 12] 



rl)'(x) ■= e - ia{x) ^)(x) 



ij'(x) 



ifj(x) 



5 a ip(x) = —iot(x)^(x) 
5 a ip(x) = ia(x)^(x) 



(3.3) 



A'^x) := A^(x) +d fl a(x) 



df,a(x), 
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where a{x) is a real function, and d~ a an operator which acts on fields and function- 
als by performing an infinitesimal gauge transformation. The condition of gauge 
invariance of our model is 

5 a e W.^]) = 

without the integral, i.e., the path integral measure times the weight should stay 
invariant. In the case of the U(l) gauge transformation ()3.3jl the measure remains 
invariant, therefore the gauge invariance reduces to a condition on the action, 

S a S[^,A] =0. 

Consider first the gauge variation of the fermionic action, 5 a S[ip,iJj}. If a(x) = const 
this transformation probes global charge conservation: this is satisfied as long as 
each term in the action has the same number of fermion annihilation and creation 
operators. But if we admit an a(x) varying in space-time, the action will in general 
change proportionally to (9 M a(x), 

5 a S[4 ! ,^j} = (f,d,a) + 0(a 2 ), (3.4) 

where we have called the coefficient of the d^a term j^. This is the motivation to 
introduce the gauge field A^. if it is coupled to j M by the gauge variation 

of the A^ term will cancel the gauge variation of the fermionic action (|3.4|) to first 
order in a. 

This raises the question how to construct the current operator j M in the action 
()3.2|) . In the continuum this is achieved by the procedure of minimal coupling. On 
the lattice, one has to use lattice gauge theory which is technically quite different 
(for an introduction, see for instance |Zinn- Justin 20021 chapter 34]). Note that (i) 
the general arguments about symmetries and Ward identities can be seen already 
in the technically simpler continuum model, and {it) the lattice formalism will be 
used in this work only in the special case of a Id lattice where the field-strength 
tensor vanishes. Therefore, I derive the Ward identities pedagogically only for the 
continuum case but give the final results also specifically for the Id lattice (in section 
14. 4|) after the lattice formulation has been introduced in Chapter EJ 

3.1.2 Minimal coupling in the continuum 

One way to construct a gauge-invariant continuum action is by minimally coupling 
the gauge field to the fermions: in the quadratic part of the action we replace the 
partial derivative by the covariant derivative 

Dy. := d^ + iA^x), 



29 



3 Ward identities in the functional RG 



such that 

£>; = e- ia D^e ia , 5 a D, = id,a. 

Likewise, the inverse bare propagator Q(d fJj ) is replaced by Q[A] = Q(D^) which is 
also gauge covariant, 

QPU = e-* a Q(D,)e ia , (3.5) 

such that the kinetic term is gauge invariant, 

Q{D'^') = Q(D^) , S a Q{D,)i,) = 0. 

Then, the current defined as the A dependence of the kinetic term, 

makes the quadratic part of the action gauge invariant. Expanding both sides of 
equation (J3.5)) to first order in a, 

where the commutator a(x)Q(x, y) — Q(x,y)a(y) = [oc(x) — a(y)]Q(x,y), and we 
have integrated by parts. 1 

Consider now the interaction term. The density-density interaction is explicitly 
gauge invariant, 



5 a V [ip,ip] = S <*J dx dyV(x - y)ip(x)ip(y)ip(y)ilj(x) 

— i J dx dyV(x — y)ip{x)ip{y)ip{y)ip{x) [a{x] + a(y) — a{y) — a(x)} 
= 0. 

On the other hand, for a bare interaction which itself is not gauge invariant (S Q Vo ^ 
0) some further counter term (j v , A M ) := — 5 a Vo is needed to make the action gauge 
invariant. Combining both the kinetic and the interaction parts of the current op- 
erator, 

f{x) :=jg(x)+j v (x), 

the action ()3.2jl will be gauge invariant. 



1 In this chapter we assume boundary conditions such that the integration by parts yields no 
boundary term. 
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Example for quadratic dispersion 

Let us give a specific example: for nonrelativistic fermions of mass m in the con- 
tinuum, the dispersion relation is linear in time and quadratic in space, Q(c^) = 
-d t + V 2 /2m. We assume a gauge-invariant interaction V$. Minimal coupling to 
the external electromagnetic potential A^x) = (<p, —A) is given in real space by 
dt d t + <p{x), — iV i— > — zV — A(x), hence the current which makes the action 
gauge invariant is constructed as 

S[4>,$,A] = J dx$(x)[-d t -<p-rt(-iV -Af]^{x)-V^^} (3.7) 

= (^,Q^)-(f,A,)-V [^^] 

paramagnetic diamagn. 

with the density operator n(x) := ip(x)ip(x). The current j has a part independent 
of A which is called the paramagnetic current and a part proportional to A which is 
called the diamagnetic current. 



3.1.3 Functional derivation of Ward identities 



The current operator j 11 [if), if), A] constructed in this way is coupled to an external 
field A^x). Thereby, all generating functionals depend on A as a parameter, 

Z[r),fj,A] := J dii Q ty,Tp) e ~^' A ,)-Vom-^,v)~(v,f) 

with measure (|2.9|) . The source fields t], fj transform under gauge transformations as 
S a r](x) := — ia(x)r)(x), S a fj(x) := ia(x)fj(x). (3.8) 

Then, the gauge variation of Z to first order in a, 



5 a Z[r],r],A] = 1^,77, Aj + ^s a rj, 



S_ 
6fj 



-tar), — I + I tar], — I + d M a, 



Sr] 



5r) 



SA, 



Z[ij,ij,A] 
Z[r),fj,A] 



0. 



(3.9) 



vanishes because Z is gauge invariant by construction. Substituting Z = e 6 we 
obtain the functional Ward identity for the connected Green functions, 



a,d. 



sg 

SiA, 



sg \ , /_ 5g 



(3.10) 
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By the Legendre transformation (|2.11j) in the presence of the external field, 

ry>, 4>, a] + (0, Q[A]cf>) := g[r), fj, a] + (0, v ) - (fj, 0) , (3.11) 

we rewrite equation in terms of the Ipi functional r[0, 0, A], 



«,^^-{r+(0,Q[A]0)} 



.a\^ + Q[A]<f> 



5T 



^,a<^ + + (0, [a,Q[A]]0) . 

By equation ()3.6|) we obtain the Ipi Ward identity 



^SiA, 



6T \ (- 5T 

— , a0 + 0, a— j 



(3.12) 



By amputating bare propagators C (without external field) from the external legs 
of equation (j3.10|) using (j2.14j) we obtain 



a, 



SV 

5iA. 



c ir-x 
°x 



uQx + Q X, a 



CiT + X 
(>X 



and by rearranging terms the Polchinski Ward identity 



a 



(x, [<x,Q]x) 



(3.13) 



On the RHS of equation ()3.13|) the non-amputated external legs are shifted. 



3.1.4 Momentum- and real-space formulation 

All of the above Ward identities still contain an inner product with the arbitrary 
gauge-transformation parameter a- q . In Fourier space the transformation ()3.8j) reads 



Vk = Vk~iY^ a -gVk+g + 0{a 2 ) 



Vk = Vk 



(3.14) 



where the mode ot- q subtracts momentum q from the fermion. The coefficient of 
a_ 9 in the Ward identity for the connected Green functions, equation (|3.1U|) . is for 
example 



ESQ 5Q 
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where the indices k± denote k ± q/2. The Ward identities for the other functionals 
has an analogous form. Taking only the one-particle component of this functional 
Ward identity by applying 1^ = ^ =0 we obtain the Ward identity given as an 



example in equation (j3.1j) . Generally, the Ward identities can be represented dia- 
grammatically as 



/ 




E 





Q 



AG! 



(3.15) 

where the wiggly line on the lhs denotes the response function of which the di- 
vergence is taken. As momentum q is transferred to the external field, momentum 
conservation implies k\ + ■ ■ • + k m = q + k[ + ■ ■ • + k' m . The Green functions in the 
large parentheses are without external field, so any one of the external legs has to be 
momentum-shifted to satisfy momentum conservation. The rightmost diagram with 
the dotted external legs shall abbreviate all the combinations of momentum-shifted 
external lines in the middle. 

Generally, for the case of global gauge invariance there is only the q = mode of 



OL^ q and the Ward identity 



- 6 1_ 

5f] 



(3.16) 



implies only that there are as many creation as annihilation operators in every mono- 
mial in the action, independent of their momenta or positions. 

In real space, the coefficient of a(x) in equation ()3.1()j) fulfills the Ward identity 



0,, 



6Q 



SiA^(x) 



JQ_ 

5r](x 



■Tj(x) + T](x] 



sg 

8fj{x) 



3.2 Cutoff Ward identities 

Generally, a momentum cutoff breaks local gauge invariance. The cutoff divides the 
fields ip into high and low modes. However, a local gauge transformation a(x) which 
is multiplicative in real space, 5 a ip{x) ~ a{x)ip{x) from equation (J3.3)) . becomes a 
convolution in momentum space, S a ipk ~ J2 q ^-q^k+q from equation (J3.14|) . This 
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shift of the momenta spoils the division of modes and poses a problem when treating 
gauge theories using flow equations. These problems would be solved if we could 
somehow define a gauge- invariant cutoff propagator Q A [A|. However, in gauge theo- 
ries, where one integrates over fluctuations of the gauge field A in the path integral, 
this is difficult except for pure gauge theories |D'AttanasiofcMorris 1996] . Note that 
these problems are absent for global gauge invariance (no momentum shift) and for 
the temperature and interaction flow schemes which have no momentum cutoff. 

There are two ways in the high-energy physics literature to deal with this problem: 
(i) One can give up gauge invariance at intermediate cutoff scales but try to ensure 
that the full Ward identities are recovered as the cutoff is finally removed (A — > 
0); this approach of modified Ward identities is explained in section 13.2.11 (ii) 
One can manifestly satisfy Ward identities in the presence of a cutoff with the help 
of an auxiliary field, the external background gauge field A (section 13.2.2)) . Both 
approaches become cumbersome if one has to truncate the flow-equation hierarchy. 
However, if one uses a manifestly gauge-invariant formulation (section I3.2.3j) such 
as the temperature-flow scheme, we show that the Ward identities can be satisfied 
exactly even in truncated flows (section |3""2""""J) . 



3.2.1 Modified Ward identities 

The Ward identities ()3.10|) which hold without cutoff are broken at intermedi- 
ate cutoff scales (A > 0) by modification terms, leading to modified Ward iden- 
tities (mwi). These modification terms vanish in the limit A — > 0. Compatibil- 
ity of flow and mwi ensures that the full hierarchy of flow equations satisfies the 
mwi at all scales A. This is sufficient for proving perturbative renormalizability 
of QED; the mwi for the connected amputated Green functions were derived in 
|KellerfcKopper 1991} KellerfcKopper 1996| and used to prove bounds for the full 



hierarchy of Green functions (see also |bdm 199*1] . and |Ellwanger 1994| for mwi in 

the Ipi scheme). For practical computations, however, the flow equations need to 
be truncated, and the mwi (and ultimately the original Ward identities at A = 0) 
are satisfied only to truncation order. Alternatively, if there are only a few relevant 
components of the flowing vertices, one can determine some of them not by the flow 
but by the mwi at every scale A, thus satisfying the mwi exactly |ehw 1 996 . 



Modified Ward identities for connected Green functions 

The definition of the cutoff connected Green functions ()2.23)1 with the addition of 
the external field can be split into a gauge- invariant part (without cutoff) and a part 
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containing the cutoff function, 

-g A { v ,rj,A} = f U,Q\AW)-V W,$]-($,r,)-(firt) det Q (j>,\Q A -QU>) (o 17% 

J detQ detQ A V ' ' 

Performing the gauge transformation (|3.3j) . (|3.8j) on all fields, the first part remains 
invariant while the latter gives to first order in a a term involving [Q A — Q] which 
breaks gauge invariance at A > 0, 



$ ae -9 A iv',v',Ai] = / l»WJ c {^,Q\A\^)-v \i>M-(i>,rt)-{n,i>) (3.18) 
J det Q A 

x ($,i[ a ,Q*-Q]i/>) e &iQ A -w) 

= ^(^^f]e~ gK ^. (3.19) 
\5r] SrjJ 

The 77, fj derivatives connect two legs of e^^ A with a propagator 
X A := [a, Q A - Q] 

which vanishes as A — > or for homogeneous a(x) = a. On the other hand, expand- 
ing Q A [r]',fj', A'] in the arguments to first order in a as in equation (|3.9jl yields 

S a g A [ V ',fj',A'} = U-ia V , j-\ + (m-^j + (^^j^^rl, A]. 

(3.20) 

We combine ()3.19j) and (J3.2()j) into the modified Ward identity 



o &G A \ T fyA^ A V/« A ^£ A \ f SQ A \ (_ 5G A 



(3.21) 

which agrees with the usual Ward identity l|P7Tjl in the limit A -> (X A -> 0). The 
modification terms have the same structure as the RHS of the flow equation ()2.24j) 
with the replacement Q A t— > X A . 

Modified Ward identities for 1PI vertex functions 

Using the Legendre transformation with cutoff (|2.27|) we can rewrite equation (J3.21|) 

as 



«, d»^J - Tr (X A [Q A + - 1 ) + (0, X\ 
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and finally obtain the Ipi modified Ward identity 




The modification term is due to the inverse bare propagator Q A —Q+Q[A] in equation 
(13.17)) . where the cutoff acts only on the A-independent part, while the A dependence 
is treated as part of the interaction (self energy), Q A — £ A [A|. The modified Ward 
identities for the Ipi vertex functions T A may be written diagrammatically as 

G A X A G A 




A7 A 

where the dashed lines denote the propagator G A X A G A and the dotted lines denote 
a momentum shift on any external leg as in equation (|3.15j) . The wiggly line in 
the third diagram means that the expansion of the inverse second derivative yields 
further one-loop diagrams with an appropriate number of vertices. 

The complete lhs, i.e., the response function together with the modification terms, 
will be interpreted in section l3.2.3l as the divergence of the effective response on scale 
A. In terms of this new response function the Ward identities will be satisfied without 
modification even after truncations (section 13.2.4)1 . 



Modified Ward identities for connected amputated Green functions 

By amputating C A from the external legs of Q A in equation ()3.21)) we obtain the 
Polchinski modified Ward identity 



36 



3.2 Cutoff Ward identities 



— , [C A aQ + D A Qa]x) + (x, [QaC A + *QD A } — ) 



3.2.2 Background-field method 



The background-field method for qed was developed in |ReuterfeWetterich 1994] 



and summarized in |FreirefeWetterich 1996] : the primary goal is manifest gauge 
invariance of the Ipi vertex functions at every scale A such that one only has to 
consider the relevant gauge-invariant couplings and not the much more numerous 
counter terms breaking gauge invariance. This comes, however, at the price of in- 
troducing an auxiliary external background gauge field A in addition to the internal 
fluctuating gauge field A. The gauge transformation of A is designed to cancel the 
modification terms in the mwi, such that T A [(p,(j), A, A] is gauge invariant under 
simultaneous gauge transformations of all fields. 

Still, gauge invariance in all fields is not sufficient to guarantee gauge invariance 
in the physically relevant ip, ip, A fields: this has to be required separately by the 
background-field identity constraining the A dependence of the vertex functions. 
Moreover, the gauge-fixing term acquires a complicated scale dependence. If there 
is no fluctuating A field, however, there is no gauge-fixing term, and in special cases 
a simple gauge-invariant construction is possible. 

3.2.3 Manifest gauge invariance 

The discussion of the background-field method (cf. section 13.2.2]) raises the question 
whether we can construct a gauge- invariant cutoff kinetic term Q A [A] in the pres- 
ence of an external field A. For a Lorentz-invariant model, —iD^ is a hermitean 
operator and so is Q[A] = Q{D l f)^ with real energy eigenvalues e(A) depending on 
the particular configuration of A. This allows to define a cutoff x A (|e(^4)|) i n terms 
of this energy, instead of frequency or momentum, by 

Q A [A] := - Q[A] =- (3.22) 

for normal operators Q[A]. For finite matrices and discrete spectra, Q[A] can be 
diagonalized and the cutoff applied separately for each eigenvalue e^, Q A (ek) : = 

Because Q A [^4] is gauge invariant by construction, equation ()3.6|) now holds also 
with cutoff, 

(«>^ 6 -^r) = MV]]- (3-23) 
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In the background-field method the background field A is coupled via the term 
(Q A [^4] — Q[A}) in the action. Using equation (J3.23J) . the gauge transformation of 
this term, 



S a A^ -^{Q^A} - Q[A}} 
oA^ 



= -iX f 

A=0 



indeed cancels the modification term X A in equation (|3.19jl . Hence, the A-dependent 
cutoff propagator restores the Ward identities without modification. 

However, in this work we concentrate on non-relativistic applications with a disper- 
sion relation as in equation ()3.7|) . which contains a real time derivative or imaginary 
frequency. Then the situation is different: Q[A] is not a hermitean operator any 
more, and it is not even a normal operator if the commutator 



[Q*[A], Q[A]} = -2dt (<p + ^(*V • A + iA ■ V + A 2 )) 



does not vanish, i.e., if the electromagnetic potential A depends on time. Then 
Q[A] and are not diagonal in the same basis, and the above definition ()3.22j) 

of Q A [A] is not applicable. Unfortunately, this condition excludes the important 
case of a finite frequency transfer u to the external field, even in the limit uj — > 0, 
which is essential for transport. Thus, energy- momentum cutoffs appear not to be 
useful in constructing a gauge- invariant Q A [^4] for non-relativistic models at finite 
temperature. If Q couples only to the homogeneous q = mode of A, the momentum 
transfer is zero and Q[A] remains diagonal in momentum space, such that one can 
use a momentum cutoff, Q A [A](K) := Q[A](K)/x A (K), for all frequency shifts ui. 

Consider, therefore, alternative flow schemes without an energy-momentum cutoff. 
If the model is regularized by finite temperature, the temperature- and interaction- 
flow schemes (cf. section l2.2.1j) allow a trivial definition of Q A [/4]. For example, in 
the interaction-flow scheme with g — ... 1 we can define Q 9 [A] :— Q [A] / g for any 
momentum transfer to the external field. It is essential that the current coupling to 
the external field is rescaled by the temperature or interaction strength just like the 
A-independent quadratic part. 

In those cases where a gauge-invariant construction of Q A [A] is possible we can 
define a new gauge-invariant generating functional 2 as 



e 



detQ A 1 ' ; 



2 Equations (|3.24|) and l|3.26|l are equivalent to equation (4.10) in [RcuterfcWcttcrich 1994 without 
fluctuations of the gauge field, a = 0. 
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which differs from the previous definition ()3.17j) in that the cutoff acts also on the 
A-dependent quadratic part. This yields the Ward identity 



a, d. 



sg 



gi,A 



SiA, 



sg 



gi,A 



Si] 



OCT] 



77, a- 



sg 



gi,A 



Sr] 



(3.25) 



on any scale A. Instead of the Legendre transformation (|2.27|) we can now define 

r gi ' A [0, 0, A] + (0, Q A L4]0) = ^[77, /7, A] + (0, r/) - (77, 0) . (3.26) 
Then r gl,A evolves under the flow equation 



<9 A r gi - A = Tr Q A C A - Tr Q A [A] (q a [A) + 



and satisfies the Ward identity 



ST^ A 



ST gi,A \ / 6T gi,A 

———,a(p + 0, a — -j— 



(3.27) 



(3.28) 



This new functional converges to the original one, r gl ' A — > T in the limit A — > 0, and 
pgi,A _ y a f or ^4 — 0. However, it has the advantage that its vertex functions are 
manifestly gauge invariant during the whole flow, i.e., the Ward identities are not 
modified. Even for a truncated flow-equation hierarchy this remains true as I will 
show in section 13.2. 41 while the flowing response vertex at A = approximates the 
exact response vertex only to truncation order, it still satisfies the Ward identities 
exactly to all orders. 



Gauge-invariant response vertex 

Let us give an explicit example of the current-response operator with cutoff. We 
diagonalize apply the cutoff on the eigenvalues and take the A derivative to 

obtain the current. To leading order in A, all modes A(q) couple independently, 
and it is sufficient to consider a single q mode. We therefore capture the generic 
situation by considering a two-state system with momenta k ± q/2 = ± and an 
inverse propagator 

- (_«♦,. + o(4 

which is diagonal for A = 0, and where the coupling to A transfers momentum q 
between the two states, J A = J tl (q\ +; -)A M (g), (J A)* = J^(— q; -; +)A IM (—q), and 
the current 

w= ( J"(g;+;-)\ 

-^(-g;-;+) )' 



39 



3 Ward identities in the functional RG 



Generally, Q[A] does not need to be a hermitean matrix but only a normal matrix; 
then the general formula (|3.22|) has to be used. For a pedagogical derivation, however, 
we assume that Q[A] is hermitean, as for instance with a purely quadratic dispersion 
Q{kfi) — \k\ 2 /2m. Then Q[A] is diagonalized, 




JA 



= ( Ly Q+ : Q - ){T X ) ( (ja)* +■ 



i 



where notably, the eigenvalues depend only quadratically on A. Replacing the diag- 
onal Q by Q A = ^ Q Q + Q ° A J , we obtain 



-JA 



and the cutoff current operator reads 

~Q + -Q- ' 

i.e., the original current rescaled by a number. As a check, the current operator thus 
constructed has the right limit, J M ' A — > J M as A — * 0, and satisfies g M J p ' A = — Q A 
as a consequence of the gauge invariance of Q A [A], equation (|3.23j) . 

3.2.4 Ward identities in truncated flows 

If the model and flow scheme permit the construction of a gauge-invariant cutoff bare 
action as in section 13.2.31 i.e., if <3 A [^4] an d the bare interaction are manifestly gauge 
invariant at any scale, then the full flow-equation hierarchy satisfies the unmodified 
Ward identities on all scales. We shall see that in this case even truncated flows can 
satisfy the Ward identities exactly. 

The most commonly used truncation of the flow-equation hierarchy without ex- 
ternal field is to set the flow of higher Green functions to zero, 

<9 A G A (A = 0) := Vm > m , (3.29) 

for some m > usually determined by practical considerations and justified pertur- 
batively in the renormalized interaction. I will now show that the Ward identities 
are satisfied if we demand the same truncation (|3.29j) also with external field, 



<9 A G A (A) := Vm > m and \/A (3.30) 
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and take derivatives with respect to in order to obtain the truncated flow equa- 
tions of the response functions G^[{A) := 8a G^(A). In particular, also the re- 
sponse functions with m > m do not flow, 

d A G^{A) := Vm>m . (3.31) 

The response-function Ward identities (J3.25)) are of the form 

d,G^(A)=zSGi(A) Vm (3.32) 

homogeneous in the number of fermion lines, where S represents the momentum 
shift on the external legs. By the truncation, neither side of equation ()3.32|) flows 
for m > mo, so the higher Green functions remain at their initial condition — given 
by the bare action — which by construction is gauge invariant and satisfies the Ward 
identities. 

For m < m , the RHS of the truncated flow equation is built up completely from 
bare propagators Q A [A] and single-scale propagators Q A [A] which by construction 
are manifestly gauge covariant on any scale A, as well as from Green functions 
G^(A) which we assume to be gauge covariant on a particular scale A. Under a 
gauge transformation, all Green functions and propagators acquire phase factors 
which cancel on all internal lines, leaving only the phase factors on the external legs 
which, in turn, imply gauge covariance of the RHS. Then, by infinitesimal induction 
the Green functions G^ dA (A) at an infinitesimally lower scale A — dA will also be 
gauge covariant, and hence the Ward identities are satisfied during the complete f RG 
flow. 

This argument is valid in any f RG scheme where the Ward identities are homoge- 
nous in the number of fermion lines, which holds in the Polchinski, Ipi and Wick- 
ordered schemes presented in this work. 



3.3 The role of self-consistency 

The generating functional Q ShA [q,f], A] defined in equation ()3.24jl has a redundant 
parametrization: there are two different ways to obtain the response functions, either 
by taking a derivative with respect to A^ or by inserting S< ^^ ipj into the path 
integral. Assuming from now on a gauge-invariant interaction Vq and dropping the 
label "gi" , we have 

J_ P -<3 K [ri,fi,A] _ f Wj] f j 8Q A [A] \ ($,QA[AW)-VoW,®-($,v)-(ri,il>) 
SA, 6 J detQA V SA» V J 

_ (5_ 5Q^4] 6\ Ah . A] 
V<V 5A, 5f}J 
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The second derivative yields terms linear and quadratic in Q 1 



5Q A 




( SQ A 


5Q A [A] 5G A \ 


SA„ 


V SAp SrjSfjJ 




5A^ Sfj ) 



(3.33) 



This self-consistency equation expresses the response functions in terms of higher 
Green functions with a loop closed, and tree terms. Equation ()3.34j) illustrates the 
relation between one- and two-particle connected Green function diagrammatically, 





+ 



(3.34) 



q,Gt: A G A G A G A 

where the dashed lines are 5Q A [A]/SA^. This relation is generally broken by trun- 
cated flows, as we shall see in section 13.3.11 Using the gauge transformation of the 
quadratic part ()3.23|) we can write the divergence of equation ()3.33|) as 

5G A \ m / r ^,^s 2 g A 



a,d L 



5iA„ 



Tr([a,QV]] 



5r] 5r] 



5rj 



6r) 



(3.35) 



and replace the response function in the Ward identity (|3.25j) to obtain the self- 
consistent Ward identity 



Tr([a,Q A [A}] 



5 2 g A 

Sr) 5fj 



Sr] 



5r] 



5g A 

5r] 



arj 



rj, ai- 



SO 1 



Srj 



(3.36) 

Note that we would arrive at the same self-consistent Ward identities to leading order 
in A if we would start with the generating functional (j3.17j) and the modified Ward 
identities (j3.21|) . By Legendre transformation ()3.26|) we obtain the self-consistent 
Ward identity for the Ipi vertex functions, 



Tr ([a,Q A [A]} [Q A [A] + 5 2 T A ) 



sr 



<f), a- 



sr 



(3.37) 



The topological structure of the A derivative on the lhs is the same as that of the 
A derivative in equations ()2.24|) . (J2.28)) . The self-consistent Ward identities for the 
Ipi vertex functions are represented diagrammatically as 



o 

Irn 
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where the dashed lines feature the propagator G A [a, Q A [A]\ G A . 

By amputating bare propagators C A [A] from the external legs in equation ()3.36j) 
we obtain the self- consistent Polchinski Ward identity 



_ (v ^ Ar ^ 2 V A \ f5V A , ^ AfA11 5V A \ (5V A \ ( 5V A 



(3.38) 

While the previous response-function Ward identities were homogeneous in the 
number of external electron legs, the self-consistent Ward identities close a loop 
on higher Green or vertex functions and thereby decrease the number of external 
legs. This "inhomogeneity" leads to severe problems: when the flow is truncated the 
lower Green functions are flowing but the higher ones are not, hence self-consistency 
is violated, and the above self-consistent Ward identities ()3.3(ij) . ()3.37j) . and ()3.38j) 
are only satisfied to truncation order. Despite intensive search a general solution to 
the problem of self-consistency remains elusive. 

In the upcoming section 13.3.11 1 will illustrate how the self-consistent Ward iden- 
tities are broken by the standard truncation that we have used in section 13.2.41 In 
section 13.3.21 1 review the conserving approximations of Baym and Kadanoff which 
are self-consistent and satisfy these Ward identities, but turn out to be in general 
incompatible with truncated flow equations except for special cases (section 13.3.31) . 

3.3.1 Self-consistent Ward identities in truncated flows 

Let us give an example of how a simple truncation breaks the self-consistent Ward 
identities. This example shows how the problem arises generically. Consider the 
Polchinski scheme with the truncated flow equations ()2.39j) . 









where slashed lines denote C A . The corresponding Ward identities ()3.38j) are 

(3.39) 



+ 



0, (3.40) 
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where dots on the external legs denote a momentum shift on any one of the legs as 
in equation (|3.15|) . and dashed lines denote [a,C A ]. This truncation is simpler than 
the ones used in practice; however, it is useful pedagogically because it demonstrates 
the problems already at second order in the renormalized interaction. A truncation 
at a higher level would display the same type of problem, only further down in the 
hierarchy where the diagrammatics is more tedious. 

If the Ward identity is to be compatible with the truncated flow, the A derivatives 
of the lhs and RHS of the Ward identity ()3.39|) for the one-particle function should 
agree: 



<9aLHs(wi) 




<9a rhs(wi) 



The second line can be rewritten using both Ward identities above and <9a [a, C A ] 

[<*,&], 



<9 a rhs(wi) 



± ± 




The lhs and RHS differ, thus we cannot complete a proof of the Ward identities by 
induction: even if the Ward identities are satisfied at some scale A, the flow violates 
them by the terms highlighted in the box. There is no reason why these terms should 
in general vanish, or why the Ward identities should be miraculously satisfied at the 
end of the flow even though they are violated during the flow. Similarly, in the Ipi 
scheme with a truncation such that the two-particle vertex function does not flow, 
the violation is of the form 







n 


-7^ 


0> 





(3.41) 



If one had not truncated the flow and Ward identity of the two-particle Green func- 
tion, they would have generated the missing terms. The difference to the formulation 
of the Ward identities in terms of response functions can be seen in ()3.32|) and ()3.39|) : 
while the one-particle response function follows the same flow equation ()3.3())1 as the 
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one-particle Green function, only with an additional A derivative, the two-particle 
Green function follows by truncation a different flow equation and is, therefore, 
not determined by the same approximation (in the sense explained below) as the 
one-particle function. Hence, self-consistency and with it the self-consistent Ward 
identities are violated. This will be illustrated more clearly in the next section 21 
on conserving approximations. 

Note, however, that this violation may be not so bad numerically: the violation 
terms highlighted in the boxes all have the structure of the terms neglected by the 
truncation, 

7^^- (Polchinski), ^CC^X^ ( 1PI ) ( 3 ' 42 ) 

but with another dashed line closed, and vice versa with the slashed and dashed lines 
interchanged. Therefore, if the truncation is justified because the truncated terms 
are small, then also the violation will be small. 




3.3.2 Conserving approximations of Baym and Kadanoff 

Baym and Kadanoff |BaymfeKadanoff 196H |Baym 1962| introduced a formalism to 
obtain conserving approximations which by construction satisfy the self-consistent 
Ward identities. We shall compare the fRG to these approximations and see whether, 
under certain conditions, the fRG might also provide conserving approximations. 



Proof of number conservation 

In order to make the reader familiar with the derivation of conserving approxima- 
tions, we shall discuss the proof of the number conservation law due to |Baym 1962| . 
During the proof it will become clear what the requirements are, and we will try 
to answer the question whether truncated fRG schemes might satisfy these require- 
ments. 

Baym and Kadanoff formulate their theory in terms of propagator lines as basic 
objects, not the source fields r], fj we have used. Therefore, they add an external 
field to the partition function, 

Z(U) := J e S[ ^H^) 

with U = U(l, T) bilocal in space-time: 3 (faUif)) = f dl dl'^(l) U(l, 1') ^(1'). The 
labels 1, 1' are abbreviations for x\, xy with all space-time components, and the 

3 This should not easily be confused with the two-particle interaction U A in Chapter 0] 
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integral J dl = J dr\ j Q 1/3 dti. Thereby, all Green functions are functionals of the 
external field U, such as the bare propagator G (l — 1'; U) and the full propagator 
G(l, 1'; U). The linear response of G(l, 1'; U) to the external field is the two-particle 
correlation function, 



L( 12, 1'2' 



8G(1,V;U) 



u=o 



[G 2 (12, 1'2') - G(l - l')G(2 - 2')] 



5U(2',2) 
The Dyson equation 

G~\l, 1'; U) = G^(l - 1') - U(l, 1') - E(l, 1'; G(U)) 



u=o ■ 



(3.43) 



shall be satisfied exactly even by an approximate G(U) and self energy E. 

The first requirement is that E be a functional of the full propagator G(U) and 
the bare gauge-invariant density- density interaction Vq. We shall see below that this 
may be relaxed a little to include any propagator that transforms in the same way 
as G{U) under change of U. 

In order to show that the approximate L satisfies the local number conservation 
law, we choose an external disturbance U that corresponds to a gauge transformation, 



dl 



^p(l) + Va(l) • + 2^Va(l)p(l)} 



The equation for the bare propagator becomes 

with the solution 

G (l, 1'; a) = e~ ia{1) G Q (1 - 1') e ia{1,) 



(3.44) 



(3.45) 



where Gq(1 — 1') satisfies ()3.44|) with a = 0. If we assume the boundary condition 
of the external disturbance to be 



a(r, r = 0) = a(r, r = —iff) 



(3.46) 



then also the solution Gq{o) satisfies this boundary condition. The full propagator 
is the solution of the Dyson equation (|3.43jl . 

{'W, ~ ^ + i + iVa ( 1 )] a } Gf ( 1 ' ^ «) - / dl I; G(a)) G(l, 1'; a) 
= 5{1 - 1'). (3.47) 
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We will now show that if £ is a functional of G(a) (as we have assumed), also G{a] 
will transform analogously to equation ([3.45)1 . To this end, substitute 



G(l,l» 



-ta(l) 



G(l,l';a)e 



ia(l') 



(3.48) 



in equation ()3.47jl and in the functional H(G(a)). If and only if all propagators in 
£ depend on U in the same way, at each vertex Vq there will be four phase factors. 
Because particles are assumed to be conserved at each vertex (by the gauge invariance 
SoYo — 0), these factors cancel exactly at each interaction vertex and remain only at 
the external legs, such that £ transforms as 

E(l, 1'; G(a)) = e~ ia{1) S(l, 1'; G) e ia{V) . 

Then, equation ([3.47)1 becomes 

- J dl e" ia(1) S(l, 1; G) G(l, 1'; a) e ia{V) 



a) e 



ia(l') 



ia(l) 



i^- + —V\\G 

dti 2m 1 



EG 



,10(1') 



5(1-1'). 



Because of the 8 function on the RHS the two phase factors cancel, and we obtain 

which is the Dyson equation ()3.47|) defining G(l — 1'; a = 0). Because the solution 
to this equation is unique, 

(5(1,1'; a) = G(l - l';a = 0) , 

and using equation ([3.48)1 . 

G(l, 1'; a) = e~ ia{l) G(l -l';a = 0) e ia{V) , (3.49) 

which is the same transformation law as for the bare propagator, equation ([3.45)1 . 
Now we expand both sides of equation ([3.49)1 to first order in a: 



d2{^^L(12, 1'2) + Va(2) 

L ut2 

i[a(l) - a(l')]G(l - !')■ 



2im 



-L(12,r2' 



2'=2 



47 



3 Ward identities in the functional RG 



Integrating by parts on the LHS using the boundary condition (|H.4fi|) and comparing 
the coefficient of a (2) yields the number conservation law for L, which is equivalent 
to the self-consistent Ward identity ()3.36j) : 



AL(12,1<2) 



v 2 - V 



lira 



= -i[5(l-2)-5(l' -2)}G{1-1') 

2'=2 



The conservation law in the 1, 1' variables of L follows if we demand in addition 
that L be symmetric in 1, V <-> 2, 2': 

<5G(1,1') 5G(2,2>) 



SU{2',2) 6U(V,1)' 

This requirement of vanishing "curl" implies, except for pathological cases, that there 
exists a functional W(U) such that 

O (M0 6W 



su(v,iy 

The <fr functional 

For completeness, let us mention that Baym and Kadanoff express the 1, V conserva- 
tion of L also as a condition on S. S is assumed there — in contrast to this work — to 
be a functional only of G(U) and Vq but not of other propagators like Gq(U), or U 
directly. Then they derive that a similar vanishing- "curl" condition must be required 
of S, 

5E(1,V) 5£(2,2') 



6G(2',2) SG(V,iy 
and hence, there exists a functional <&[G(U), Vq], such that 
5$ 



E(1,10 



6G(V,iy 



and S'E/SG is the effective particle-hole interaction. In equilibrium when U(l, 1') = 
S(h - t[) U(r u r[), W and $ are related by W = $ - tr(SG) - trln(-G). An 
approximation for that can be written as 5&/5G is called ^-derivable. 



Diagrammatic interpretation of conserving approximations 

Let us give a diagrammatic illustration of the Baym-Kadanoff formalism. In the 
exact theory W(U) contains every vacuum diagram of perturbation theory with bare 
vertices Vq and propagators Gq(U). A derivative with respect to U acts on the Gq(U) 
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lines in each diagram by plucking out one line. The full propagator G(U) = SW/SU 
is, therefore, made up of all connected vacuum diagrams with two external legs. A 
second derivative with respect to U plucks out a second line from each diagram, 
generating L = —5G/5U, which consists of all diagrams with four external legs. 

Consider approximations defined by choosing a possibly infinite subset of all Feyn- 
man diagrams for W(U), G(U), and L. In principle, one could choose completely 
different sets of diagrams for each correlation function. But an approximation is con- 
serving if and only if L = —5G/5U is satisfied exactly, i.e., the diagrams contributing 
to L are obtained by plucking out one line in each diagram of G, or equivalently, 
the G diagrams are obtained by closing one loop on each diagram of L. Such an 
approximation for L is conserving in the 2, 2' variables; it can be made conserving 
also in the 1, 1' variables if we demand further that the diagrams of G are derived 
from those of W by plucking out one line. 

Obviously this leaves no room for truncations only on a certain level of the hier- 
archy of Green functions: the flow equation for the zero-particle component W(U) 
must completely determine all higher flow equations in order that the classes of 
diagrams contributing to each Green function are compatible. 

3.3.3 How important is self-consistency? 

It has become clear that at the core of the problem are not the response-function 
Ward identities but self-consistency: the relation L = —5G/5U is not satisfied in the 
above example where the one-particle function flows but the flow of the two-particle 
function is truncated. This incompatibility of the truncated flow with conserving 
approximations raises the question under which circumstances it is a problem not to 
have a conserving approximation. The answer depends on the physical problem at 
hand: there are problems which work surprisingly well in the truncated fRG, while 
others fail miserably. 

A favorable example are the Id impurity problems presented in Chapters |U and El 
Already for simple truncations the fRG results are very close to the exact asymptotic 
solution known from Bethe ansatz. This suggests that the terms neglected by the 
truncation are small. Following the discussion surrounding equation ()3.42|) . also self- 
consistency is then fulfilled to a high degree of accuracy. Moreover, in the model 
of Chapter 0] the self-consistent Hartree-Fock approximation, a simple conserving 
approximation, is known to produce the wrong physical phase (charge-density wave). 
This is due to the fact that an approximation, while being conserving, can miss 
important contributions. Generally, if the truncated flow does not diverge, one can 
even determine some components of G(U) and L self-consistently: one can use any 
truncated flow to obtain an approximate W A (U) in the presence of U and take 
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numerical derivatives with respect to U to obtain approximate values for G(U) and 
L which are conserving by construction (cf . section 14. 3 j) . 

An example of the opposite situation where the violation of the self-consistent 
Ward identities is disastrous is the reduced BCS model |shml 2()()4j . This model is 
solved exactly by the self-consistent Hartree-Fock approximation. In the presence of 
a tiny symmetry-breaking term (gap) of magnitude e in the action, the interaction 
grows very large to 1/e but does not diverge. In the truncated flow equation, this 
result is only reproduced if one has exactly the correct value of the gap (a component 
of the one-particle function) in the flow equation for the two-particle interaction. 
This would be guaranteed by the self-consistent Ward identity relating the one- and 
two-particle functions; however, even if the truncated flow violates the Ward identity 
only slightly the interaction may diverge prematurely at A > 0. This problem can 
be solved by a modification of the flow-equation hierarchy [Katanin 2004] : if the 
single-scale propagator S A in the Ipi scheme is replaced by the A derivative of the 
full propagator —d^G K = G K {d^Q A — dfjl )G A in the truncated flow equation for 
the two-particle vertex, then the fRG flow reproduces exactly the self-consistent 
Hartree-Fock solution for models where the two-particle interaction has a reduced 
momentum dependence. This modification leads to significant improvements also in 
the single-impurity Anderson model |hmps 2004] . 

There are promising approaches to self-consistency by writing the fRG flow equa- 
tions in terms of both fermionic and bosonic degrees of freedom. For instance, (i) 
jWetterich 2002] thereby obt ains the flow of 2pi vertex functions, (ii) The Luttinger 
model without backscattering is treated in |sbk 2004] ; using the separate conserva- 
tion of the number of left and right movers the flow-equation hierarchy can be closed 
and solved exactly. 

3.4 Summary 

In this chapter I have shown how Ward identities, which express the symmetry of 
the Hamiltonian in terms of Green or vertex functions, are derived in the functional 
formalism. There are two common formulations of Ward identities. In the field- 
theoretical and high-energy physics literature, Ward identities are written in terms 
of response functions. This form of the Ward identities relates for example the 
self energy to the current response but assumes no particular relation between the 
self energy and the interaction. A momentum cutoff breaks these response-function 
Ward identities. This is a problem for the treatment of gauge theories, leading 
either to modified Ward identities or the introduction of a background gauge field 
A. Alternatively, one can use a manifestly gauge-invariant flow scheme such as the 
temperature flow. It is shown that in this case even truncated flows satisfy the 
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unmodified Ward identities exactly on all scales. 

In the condensed-matter literature, a different form of Ward identity is more com- 
mon which assumes self-consistency: the response functions in the Ward identities 
can be expressed as higher Green functions with a loop closed by a special propaga- 
tor. An example of approximations which satisfy these self-consistent Ward identities 
without cutoff are the conserving approximations of Baym and Kadanoff. However, 
I have shown that self-consistency is violated by common truncations, which neglect 
the flow of Green functions beyond a certain level in the flow-equation hierarchy. For 
special models with exact mean-field solutions, modified truncated flow equations are 
known which reproduce these solutions. 
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4 Functional RG technique in one dimension 



In this chapter I introduce the one- dimensional lattice model of a Luttinger liquid 
(cf. Chapter^ and explain in detail how the f RG is used to compute observables such 
as the effective impurity potential, the density-response vertex and the conductance. 
This chapter is organized as follows: in the first section l4~Tl the microscopic lattice 
model is defined. In section [Ql I show how the flow equations are set up and solved, 
with an emphasis on the finite-temperature flow and a new efficient algorithm in Id. 
I proceed to explain how to compute the conductance in the f RG framework in section 
14.41 giving an argument why vertex corrections play no role in our approximation. 



4.1 Microscopic model 

Consider a model of spinless fermions on a Id lattice with nearest-neighbor in- 
teraction and various types of impurity potentials. Following [aemmss 2004] . the 
Hamiltonian has the form 

H = H + H iQt + H imp 

with the kinetic term given by nearest-neighbor hopping with an amplitude —t and 
chemical potential \x (I will henceforth choose units such that t — 1, and the lattice 
spacing a = 1), 

Ho = ~tJ2 ( c ]+i c j + c k+i) - ^ XX"' 

3 3 

where rij = CjCj is the local density operator. The nearest-neighbor interaction of 
strength Ujj + i = Uj+ij on the bond between sites j and j + 1 enters as 

H int = X U 3,3+l n 3 n 3+l 
3 

while the static impurity potential V 3 j< is represented by a term 

#im P = X V jr r j ( 'y 

3,3' 

For the conductance calculations, we will couple an interacting system on lattice 
sites 1, . . . , N on both sides to semi-infinite, non- interacting leads which are described 
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left lead 



Jl 



sample 



Jr 



right lead 



• • • 







t, 



1 2 



N-1 N 



t 



R 



• • • 



N+1 



Figure 4.1: The microscopic setup of the system with leads. The current operators Jl,r 
will only be needed for the conductance calculations. 

completely by H ffigure l4~T|) . In order that electrons do not scatter off the beginning 
of the interacting region, the interaction has to be switched on smoothly. Explicitly, 
we choose a spatial profile 



_ LT arctan [(j - j s )/w] - arctan [(1 - j s )/w] 
3 ' 3+ arctan [(N/2 — j 8 )/w] — arctan [(1 — j 8 )/w] 



(j = l,...,iV/2) 



for the left side of the system, and likewise for the right, where U is the bulk inter- 
action. We have chosen the parameters w = 4 and j s = 56 such that the interaction 
falls to 10% of its value over a typical distance 22 lattice sites. 1 



4.1.1 Projection method applied to the wire 

In order to treat this infinite system with leads numerically, we express it (exactly) 
by an effective Hamiltonian on the iV-site interacting region via the projection tech- 
nique. Consider splitting the Hilbert space of the Hamiltonian H into disjoint sub- 
spaces with projection operators P + Q — 1: 



H 



Hp P H P Q 
Hqp Hqq 



The one-particle Green function is the resolvent 
1 



G(z) :-- 



Gpp( 
Gqp( 



Gpq(z) 
g qq( z ) 



z-H 

with components |ptvf equation (2.7.23)] 

„ 1 „ 1 



Gpp(z) 



H 



:P 



zP-H 



pp 



H 



PQ 



ZQ - Hqq 



-H, 



QP 



(4.2) 



1 Although the difference \U- 



- U 



js ,js + l\ 



= :S P p( 2 ) 

0.08 U is rather large, the corresponding back- 



scattering component of the effective potential S is typically below 10 4 , and the conductance 
at T = deviates from the unitary limit by less than 10 -8 . 
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Gpq{z) = P^-frQ = -G PP (z)H PQ 1 

Z — tl Z(4 — riQQ 

G Q P(Z) = Q^-JtP = - n l rr H Q PGPP{Z) 
Z — H ZLl — Hqq 

Gqq(z) = Q^-rjQ = „ l — + ^ l — H QP G PP (z)H PQ 



ZQ - Hqq ZQ - H QQ ZQ - Hqq 

For the wire, P shall project onto the Hilbert space of the states in the interacting 
sample with site indices 1, . . . , N, while Q shall denote the remaining non-interacting 
leads on sites . . . , and N + 1, . . . . Thus, 



N-l N 



imp 



H PP = ~t^2(\n+ l)(n\ + \n)(n+ 1|) -^^|n)(n| +H mt + Hi 

n=l n=l 

H PQ = -t L \l)(0\-t R \N)(N+l\ 
H QP = -t L \0)(l\-t R \N+l){N\ 

Hqq = —t (\n + l)(n\ + \n)(n+l\) - fi ^ \n)(n\. 

n<0, n<0, 
n>N+l n>N+l 

The left and right leads in Hqq do not couple directly but only through the sample. 
In order to compute Spp in equation ()4.2j) we need to know the Green function of 
the (left) lead at the interface site 0, (0\(zQ — Hqq)" 1 ^). To this end, consider a 
semi-infinite lead ranging from — oo up to some site j, and denote the Green function 
at the rightmost site j as <7l(z) := (z — H )j-. Consider adding one more site j + 1 
to the right, with the same hopping amplitude —t. Because the lead is semi-infinite 
and homogeneous, the Green function g' L {z) at the new site j + 1 should be the same 
as gh{z). Again using the same formula ()4.2j) with Q now denoting the states on 
sites — oo, . . . , j and P for j + 1, 

[g'dz)Y l = {z- H ) j+1J+1 - (H ) j+1 j - //,;) J (#o)j,j+i 
= (z + fj)-tg L (z)t= [gdz)}- 1 . 

This leads to a quadratic equation for gi(z), 

t 2 g 2 L (z)-(z + fi) g L (z) + 1=0 

with solution (t = 1) 

g L (z) = \ [z + pi T V4-(z + /i)2) (4.3) 

and gR^z) = gh(z) because the right lead has the same structure as the left lead. 
gL,{z) has a branch cut at the real axis Imz = 0. The sign ± is chosen such that the 
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imaginary part of gi{z) changes sign across the branch cut. For the local density of 
states of the leads we obtain 

p L {u) := --lmg L (u + iO) = -^-a/4- (u + /i) 2 , 
ix An 

for | a; + fi\ < 2. Now we have all ingredients for the lead contributions to the 
propagator of the sample, equation (|4.2j) : 

T, PP (z) := Hpg— — — Hq P 

Z(4 ~ n QQ 

= t L \l)(0\(zQ-H QQ y 1 \0)(l\t L 

+ t R \N)(N + 1\{zQ - Hqq)~ x \N + l)(N\t R 
= tlg L (z)\l)(l\+t 2 R g R (z)\N)(N\ (4.4) 

G j0 (z) = (j\G PQ (z)\0) 

= -(j\G PP (z)\l) t L (0\(zQ - Hqq)- 1 ^) 

= -G, n {z)t L g L {z) (j = l,...,N) (4.5) 

G N+1J (z) = (N + l\G QP (z)\j) 

= —{N + l\(zQ- H QQ )- l \N + 1) t R (N\G PP (z)\j) 

= -g R (z)t R G Nj (z) (j = l,...,N) (4.6) 

G N+h0 (z) = (N + l\G QQ (z)\0) 

= (N + l|(zQ — Hqq)" 1 ^) (= as left and right leads do not couple directly) 
+ (N + l\(zQ - Hqq)~ 1 \N + 1) t R (N\Gp P (z)\l) t L (0\{ Z Q - Hqq)- 1 ^) 
= g R {z)t R G Nl {z)t L g L (z). (4.7) 

In the non-interacting leads, Hqq = Hq. In the sample, the quadratic (non- 
interacting) part of Hpp shall include the impurity: 

The lead contribution is also independent of the interaction, 

£ lcads (z) := £pp(z) = t 2 L g L (z) |1){1| +* R g R (z) \N)(N\, (4.8) 

and all non-interacting contributions are combined into the inverse bare propagator 
in the sample, 

QiAz) = lGo\z)]jf ■= zSj? - C,r - >-jr^- (4-9) 

The interaction H mt creates an effective one-particle potential £(z, T) at temperature 
T. The full sample propagator is determined by the Dyson equation, 

[G-\z,T)} n , = Q jf (z) - E jr (z,T). (4.10) 
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4.1.2 Bare vertices 

The f RG is formulated in terms of a path integral weighted by the action, so the bare 
and interaction parts of the Hamiltonian have to be translated into the Lagrangian 
language: the inverse bare propagator Q(z) defined in equation (J4.9)) becomes the 
quadratic part of the Lagrangian, and the bare interaction H mt becomes the inter- 
action part: 

Sty,$]=T^Q jj tiu n )$ j (iu n W j .{iu n )+ A',2';i,2^(l')^(2 , )^(2)^(l) 

n,j,j' 1,1',2,2' 

where Iv^'xi is the bare antisymmetrized interaction vertex equivalent to H [nt , 
and the sum over 1, 1', . . . includes Matsubara frequencies and lattice indices. The 
nearest-neighbor interaction conserves frequency but has no further frequency de- 
pendence. The remaining spatial dependence has the real-space form 

I j' 1 ,j' 2 ;h,h = U hJa[ S hJa-l + *iirfa+l](*JiJi<Ws ~ 5 h,j' 2 5 h,jO- ( 4 - n ) 

In the bulk Ujj + i = U is homogeneous, hence the bare interaction is translationally 
invariant and can be expressed in a momentum basis: 

h'^-MM = 2£/[cos(X - fci) - cos(^ - h)]S { ^ k2ik , i+k , 2 (4.12) 
where the Kronecker 5 implements momentum conservation (modulo 2ir). 



4.2 Functional RG flow equations 

For computing the properties of the Id fermion system, we use the Ipi version of 
the fRG, cf. section 12.2.31 We cut off the infrared part of the free propagator on 
a scale A and differentiate the generating functional for the vertex functions with 
respect to this scale. Thereby we obtain an exact hierarchy of flow equations for the 
irreducible vertex functions. 



4.2.1 Truncated IPI flow equations 

As in jAEMMSS 2UU4j . we choose a cutoff 
propagator Gq (iu) (with the notation Go 

G#(iu;) := X »GoM, 

where the characteristic function x A ( ti; ) is 
ishes on the low-energy modes. The exact 
T > 0, as explained below. 



in Matsubara frequency. The cutoff bare 
= C) is defined by 

(4.13) 

unity on the high-energy modes and van- 
form of x A depends on whether T = or 
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Second, we truncate the infinite hierarchy of flow equations by neglecting the 
flow of the three-particle and higher vertex functions. This closes the hierarchy of 
flow equations for the one- and two-particle vertex functions and is justified by a 
small renormalized interaction. The results agree remarkably well quantitatively 
with known exact results (DMRG, Bethe ansatz) not only for weak but also for 
moderate interaction strength |aemmss 2004^ [Andergassen 2005| . 

The truncated Ipi flow equations are written in terms of the one-particle vertex E 
at scale A (self energy) and the two-particle interaction vertex 2 T A |AEMMSS 2004j . 

<9 A £ A (1',1) = -T^e^ 20+ ^ A (2,2 , )r A (l , ,2 , ;l,2) (4.14) 

2,2' 

<9 A r A (l', 2';1,2) = TJ2J2 ° A ( 3 '> 3 ) S ^ 4 ') 

3,3' 4,4' 

x [r A (l / ,2 / ;3,4)r A (3 / ,4 / ;l,2) 

- r A (l', 4'; 1, 3) r A (3\ 2'; 4, 2) - (3 <-> 4, 3' 4') 

+ r A (2', 4'; 1, 3) T A (3', 1'; 4, 2) + (3 <-> 4, 3' <-> 4')j . (4.15) 

The indices 1,2,... label both frequency and spatial indices. The full propagator is 
determined from the self energy by the Dyson equation 

G A = [Q A - E A ]-\ (4.16) 

where Q A = [G A ] _1 = Q/x A is the inverse bare cutoff propagator. S A denotes the 
single-scale propagator which will select only modes with frequency near A for the 
frequency cutoff, 

S A : =G A Q A G A = -x A - l — K Q ^ —. . , (4.17) 

where the dot denotes d&. The convergence factor e tu)2 ° + in the S A flow equation is 
only necessary to define the initial condition of the flow at A = Ao — > oo (see below). 

Parametrization of the interaction 

The two-particle interaction vertex r A (l', 2'; 1, 2) is in general a very complicated 
function. However, in the low-energy limit the flow of S A and T A is dominated by 



very few channels of the interaction Andergassen 2005 , lAEMMSS 2004]. At any A, 



perturbation theory in the renormalized interaction strength provides a guide for a 
simple parametrization of T A . For a model of spinless fermions with nearest-neighbor 



2 Not to be confused with the full functional 
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interaction and only few impurities (i.e., not disordered), the following parametriza- 
tion captures the qualitative as well as the quantitative features of Luttinger liquids 
very well. The renormalized interaction vertex T A is assumed to be frequency inde- 
pendent beyond the conservation of frequency, and the RHS of the flow equation is 
evaluated with all external frequencies set to zero. Likewise, its spatial dependence 
is approximated by a renormalized nearest-neighbor interaction U A , and in the flow 
all external momenta are projected to ±fcp. The internal lines carry propagators 
without the self energy or the impurity potential, which would lead to corrections 
only at higher order in the interaction. 

In the bulk Ujj + \ = U is homogeneous, and T A has the momentum-space form 
(independent of frequency) 

n>,« 2 MM = 2f/ A [cos(A;' 1 - A*) - cos(A£ - ^C+U^ 
which is just the bulk bare interaction (|4.12jl . 

h' v k' 2 -MM = 2C/[cos(A;; - k x ) - cos(k' 2 - fci)]^+fc 2 , fei +^ 

rescaled by U A /U. The particular flow equation for U A depends on the cutoff chosen 
and will be given below for several types of cutoff. However, the general form is a 
consequence of the above parametrization, 

d A U A = -{U A f T £ *» / Y-fip, w) (4.18) 

where f(p,u) is the sum of the three different channels (PP, PH, PH') in equation 
()4.15j) jAEMMSS 2004] . For instance at half filling where fx = and £ p = -2 cos(p) = 

s— pi 



jr , = 2sin 2 (p) cos 2 (p) [1 + sin(p)] 2 

(iuj - Q(-iuj - £_ p ) (iu-Q 2 (icu - Q(-iou - Q 



— cos 2 (p) 



1 1 

+ 



(4.20) 



Going from the bulk Ujj+i = U back to the lattice, we apply the U A flow equation 
for each U A j +1 locally. 

Parametrizing T A by a renormalized nearest-neighbor interaction U A j +1 has the 
great advantage that the self energy is a tridiagonal matrix in real space: only the 
matrix elements S^.^^ are non-zero. The tridiagonal flow equations for a general 
cutoff are 



r=±l 



(4.21) 

d^ ±1 = -Tj2x A u A ±1 s A ]±1 M. 
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Note that the self energy remains independent of frequency and real because the 
interaction does not depend on frequency in our parametrization. 

4.2.2 Frequency cutoff at zero temperature 

At zero temperature we choose the sharp cutoff 

X A M := 6(M - A), xV) = S(\u\ - A) (4.22) 

which cuts off all modes with frequency smaller than A. As explained in section l2\2.31 
a sharp cutoff at T = allows to integrate over the step functions analytically 
which leaves only smooth propagators G in each loop diagram ()2.34j) . 

G A (iuj) := [Q(iu) - E A ]-\ (4.23) 

as opposed to G A (iu>) from equation (J4.16)) which has a step at \u\ = A. The flow 
equations (|4.21j) at T = for a sharp frequency cutoff then read (with \— > 



= E E U A j+r Gf +rJ+r (iu) (4.24) 

iAj =±Ar=±l 

d ^ = ^k E U A j±1 Gf !j±1 (iu) . (4.25) 



ui=±A 



Flow of the interaction 

Inserting equation (|4.22j) into (|4.18jl . the flow equation for U A is 

d A U A = (U A r±- l^-fM. (4-26) 

w=±A J 

For instance at half filling, the momentum integral over the bubble f(p,u) is 

1 1 



uj=±A 



E /(^) = -/^cos 2 (p)2Re 

dp [2 cos (p)] 4 



(iA-Q(-iA-Q (zA-Q* 



2vr(A 2 + [2cos(p)] 2 ) 2 
A 2 + 6 



1 - A- 



(A 2 + 4)3/2") 
The differential equation (J4.26j) separates, 
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Integrating each side separately from A = oo down to A with initial condition 

U A =°° = U, 

1 1 1 / 2 + A 2 
A 



and finally 

U A = 7 ^ ^77-^ ■ (4.27) 



a - m) u iw 1 - c//(2 ' r) 



Even away from half filling, the flow equation can be integrated analytically using 
contour integration |Andergassen 2005] . 

At T = 0, the simple expression ()4.27|) yields the correct low-energy asymptotics 
to second order in the renormalized vertex and moreover contains the second-order 
corrections from the lattice dispersion to the vertex (with all external lines at the 
Fermi surface) at any scale A. 

Initial conditions 

The flow, being given by the solution of an ode in A, is determined uniquely by the 
flow equation and the initial condition. At the initial upper cutoff scale A = Ao, 
the initial condition has contributions from the bare interaction and from the bare 
impurity potential jAEMMSS 2UU4j . At T = 0, only the combination G' 1 = Q - £ A 
appears in the flow, such that there is no difference whether one treats the impurity 
potential as the initial condition for X A at the beginning of the flow or, alternatively, 
as part of the bare propagator Q as in this work. The contribution of the interaction 
as A = A — > oo is 

S U' := ^E J ™ (4-28) 

2 

IWl 2 := h',2';l,2 i (4.29) 



where ii',2';i,2 is the bare antisymmetrized interaction (j4.11j) . The initial condition 
for the self energy is usually compensated by a local potential to avoid that the 
filling changes (see section 14.2.61 below). Only if the interaction extends to an open 
boundary without lead, a boundary term S A ° undary := —U/2 remains. 

4.2.3 Frequency cutoff at finite temperature 

At T > 0, a sharp cutoff x A (a; n ) = 0(|a; n | — A) analogous to the T = case would 
lead to a 5 peak on the RHS of the flow equation of S A at every fermionic Matsubara 
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frequency A = \u n \. This creates finite jumps in the integrated E , therefore, the 
G integration formula (|2.33|) cannot be applied as it is only valid for continuous 
functions. For the ode integration, a smooth right-hand side is best, but this has to 
be balanced against the number of Matsubara frequencies uj n on which x h (uj n ) > 
for every particular value of A. 




Figure 4.2: The cutoff function used at finite temperature. 
As a compromise, we use a cutoff function (cf. figure l4~! 







\Wn\ < A — ITT 



\ + ^3^ A-ttT < \uj n \ < A + vrT 
1 A + 7rT<|wJ 



with the accompanying A derivative 



A - 7rT < \u n \ < A + 7rT 







otherwise. 



A Matsubara sum over the single-scale propagator (j4.17j) with this cutoff contains 
exactly one term ±a> n in the Matsubara sum, 



\lUJ r 



The flow equation for the self energy E A has the same general form as for T = 0, 
but for T > 0, the single-scale propagator S" A cannot be further simplified to G: 



w n ~±A r=±l 



Q(iuj n ) - x A K)E A 



Q(iu; n ) 



Q(iu n ) - x A K)S A 



(4.30) 
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1 



Qiioj) 



1 



x 



Q(^ n )-x A K)s A 



Q(iu n ) - x A (^n)S 



A 



-I 3 J'±l 



At every A = cj n — 7rT, a switch from u; n to u; n _i occurs and the integrated flow has 
a kink, instead of a jump for a sharp cutoff. At large frequencies, the self energy is 
still so small that an adaptive step-size integration algorithm efficiently takes steps 
of many Matsubara frequencies at once, while at small frequencies it inefficiently 
takes many small steps to resolve the kink accurately. But since the positions of the 
kinks are known, we can instead integrate only in the interval between each pair of 
kinks. Empirically, the switch between continuous integration and integration only 
between kinks is best done around uo n pa t. As there are 0(1/T) such intervals, the 
runtime for the complete flow scales as 0(N/T). 

Flow of the interaction 

The flow of the interaction is simpler: as self-energy corrections on the internal lines 
of the interaction flow equation would be of 0((T A ) 3 ), we take only bare propagators 
on the internal lines, such that the bare single-scale propagator is (setting £ A = 0) 



which is the same as for T = 0, except that for each A one has to insert instead the 
nearest discrete u> n into the propagator. At half filling, 



At T = 0, X A and T A flow substantially only in the range O.lt < u < t. At 
T > 0, since the lowest fermionic Matsubara frequency is u = ttT, S a and T A are 
renormalized substantially only for T < 1. The interaction is still correct to second 
order in U if all external legs lie on the Fermi surface, but the frequency dependence 
at second order becomes important with inelastic scattering. 

Initial conditions 

At T > 0, the combination Q — x A ^ A appears in the single-scale propagator in 
equation (|4.3()jl . i.e., Q and S A do not enter on an equal footing. For the exact 
theory and full hierarchy of flow equations, the final result should be independent of 
the specific form of the cutoff, but in the truncated flow it does make a difference. 
Both the bare impurity potential and the lead contribution are contained in the 
bare propagator ()4.9j) , such that the self energy receives only contributions from the 
interaction, and the same initial conditions (|4.28jl and ()4.29j) as for T = apply. 



5*0 ( iuJ n) 




Q(iu n ) 
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4.2.4 Temperature flow 

In the temperature flow (see section 12.2.11 and |Honer kampfeSalmhofer 2001| ), the 
kinetic term in the action, 

n k 

is replaced after a rescaling of the fields, ip T~ 3 ^ijj, by 

iu n - i k 

n k 



U> n ,k r= 



We shall postpone the leads for a moment but include the impurity potential in £ as 
usual. With the new convention that the scalar product shall contain no further fac- 
tor of T, we define the rescaled inverse bare propagator (indicated by the superscript 
T) as 

qT = iuJ n ~ £ 

Vt 

which yields the rescaled full propagator (with rescaled self energy S T ) 
G T 1 ^ 



Q T - S T iw„ - 5 - -JTTF ' 

The T derivatives of these quantities are 

d iuj n 
dT {lUJn) = T 

d „ T iuj n wJ n — £ iu n + £ i i ;i i ■ i 

■jt q =j^~^fy^ = ^fy^ ( + Slg 

nT i i i 
s T = — =- — =— {iuj n +i) 



Note the unusual + sign in dQ T /dT which is due to the fact that one does not just 



have a multiplicative cutoff function x T = VT but the temperature also enters in 
the Matsubara frequency uj n . 

To understand the structure of the flow equation, consider first the RG flow of the 
self energy without leads, feedback of the self energy or vertex renormalization, i.e., 
first-order perturbation theory: 



S T = —C T = — 



dT dT \ ilo, 



n 
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dT 



C T 



_d_ 

dT 



J n 0+ 



d 



iuj n -ij dT \ y/T J 



The flow equation of the rescaled self energy E T is thus 



dT' 



d_ 

dT 



With the initial condition V^o 



To^oo 



tr(f//(0) 



To^oo 



\ tr(t/) as in equation 



(J4.28j) we can integrate from T = Tq down to T and obtain 
/ TS T = tr(f//(0) 



as expected from the relation between rescaled and original Green functions, equation 

Notice one peculiarity of the temperature flow: since the Green functions of the 
leads depend on the frequency icu n , they also depend on temperature and have to be 
differentiated appropriately in dQ T /dT. At half filling the expressions are particu- 
larly simple, 



Q T {iu} n ) 



leads , 



T 



d_ 

dT 



1 

2 
1 

T 



r , - i sgn(cj n ) a/4 + 
sgn(u; n 



g L [iu n ) + 2i 



d_ 

dT 



27-3/2 _ 



tLu n -^ ds (iu n )+^-At 



sgn^j 



{t 2 L \l)(l\+tl\N}(N\) 



Temperature flow of the interaction 



The temperature flow of the interaction U T has the structure of the usual one- loop di- 
agrams of the perturbation expansion in the original ip fields, but with a temperature 
derivative of the particle-hole and particle-particle bubbles |HonerkampfcSalmhofer 2001| . 
In our Id case, 



~dT~ ~ dT 



where 



(4.31) 



(4.32) 
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is the sum of the three bubble contributions. For instance at half filling, 
B T = - 



dp 2 
— cos (p) 



T 



dp |"cos(p) ^ k ( cos(p) \ 



2ix 



T 



(p) ( 



AT V 



tanh 



cos(p) 



1 



In the limit T — > 0, _B T — ► while for T — > oo, -B T — > vanishes. Equation (|4.31j) 



can be written as 



Integrating from T = oo (with f/ T 00 = U) down to T, we obtain 

if u 



1 - UB T ' 



For T = 0, the result U T = U / {l — U / (27r)) agrees with the frequency-cutoff result 
jOZD at A = 0. 



Initial conditions 

The initial conditions for E = V^E 7 and C/ = £/ T in the limit T = T 
the same as for the frequency cutoff, equations (|4.28J) and (|4.29J) . 



oo are 



4.2.5 Interaction flow 



In the interaction flow scheme (cf. section 12.2.11 or |HRAE 2()()4j ) the propagators 
are slowly switched on by a global scale factor x A — AS 9 — 0...1, irrespective of 
frequency or momentum. We define the inverse bare propagator 



Q{iuj) iuj-Y}™ ds {iu)-£ 



9 



and the full propagator 
G 9 {iu) :-- 



9 



iuo -E lcads (zu;) - f- gY.9 

where £ shall include the impurity potential V. Then the "single-scale" propagator 
(which is not at all single-scale) is 



S3 ■= G g Q 9 G 9 



Q(iuj) — gH'- 



1 



Q{iu) - gYp' 
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Furthermore, because the bare G$ = gG$ contains one g factor and S$ = Gq contains 
none, the flow equation for the interaction vertex contains the combination GqSq + 
SqGq = 2gG G in the bubble, 

^-U° = (U 9 ) 2 2gB T 
dg 

with the bubble integral B T at temperature T, equation (|4.H2jl . Integrating from 
g = (with U 9=0 = U) up to g, 

v- u 



g 2 UB T 

(T = 0). 



l-g 2 U/(2ir 
At the end of the flow, 

U 9=1 = = U T . 

1-UB T 

The initial condition for the self energy is S 9=0 = because all propagators vanish, 
but r 9=0 is the same as for the frequency cutoff, equation (|4.29jl . 



4.2.6 Initial conditions for general filling 

At any filling n, we wish to fulfill two conditions for the model without impurities: 
(i) at the end of the flow, the density profile should be uniformly nj = n, both in 
the interacting region and in the leads, such that power-law exponents depending 
on the density can be read off reliably; (ii) no backscattering should occur at the 
ends of the wire where the interaction is switched on, i.e., the transmission without 
impurity should be perfect. In order to achieve this, we allow the freedom to add 
to our microscopic model a local potential in the wire that depends only on the 
interaction strength but not on the temperature, while with changing temperature 
only the global chemical potential fi(T) may be adjusted. 

At half filling, these conditions are met without a local potential and by setting 
/j, = 0. Away from half filling, we use the following procedure: 

1. At zero temperature, the global chemical potential is defined as [t := — It cos(A;f) 7^ 
such that ^ F = 0, where the Fermi wave vector kp = nn at filling n. Second, 
finding the local potential at all iV lattice sites such that at the end of the 
flow nj = n (without impurities) is a complicated ^-parameter optimization 
problem, but we obtain a very good guess in the following way: the problem is 
first solved in the bulk where the global potential, implemented as the initial 
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condition of the self energy £ b ° lk , is tuned self-consistently until = 0. 

Then the initial condition of the self energy on the lattice (the local potential) 
is defined as 

jj '~ 2U bulk 

such that in the middle of the system where E/jj+i = U is homogeneous, the 
bulk initial condition holds, while towards the ends when the system becomes 
non- interacting Ujj+i — > 0, no local potential is necessary. Equivalently, the 
local potential may also be written as part of the interaction term (|4.1j) . 

#int = Yl U Jj+l( n 3 ~ Z/ )( n J+l - V ) 
3 

where v = (1 - S^ Ik /C/)/2. 

2. Now the temperature T > is switched on and /i is tuned self-consistently until 
the bulk density n(/i) = j- J dk f{e^ — fi) = n at the end of the flow, with the 
non-interacting bulk dispersion e& of the leads. By the above construction also 
the density in the interacting region, obtained using the dispersion e& + with 
renormalized hopping and potential, has the same value to a high accuracy, 
and the density computed using a flowing density vertex (see below) deviates 
by less than 1% for \U\ < 1. 

Only at this stage, impurities are inserted into the lattice and change the homoge- 
neous density profile. 

4.2.7 Algorithm for tridiagonal matrices 

With the parametrization of the interaction by one single bulk parameter U A , the 
most time-consuming part of the ODE flow is the flow of the self energy. On the RHS 
of the E A flow equation, both Q and E A are tridiagonal matrices, hence it involves the 
inversion of tridiagonal matrices. I have developed an efficient inversion algorithm to 
compute the RHS in O(N) time at any temperature. Therefore, the runtime scales 
only linearly with the system size N, and the self energy of systems as large as 
N = 10 7 sites (at T = 0) can be computed. At T = 0, only the tridiagonal part 
of the inverse tridiagonal matrix G A is needed in (|4.24|) . and an algorithm by the 
author was presented already in |aemmss 2UU4j and is reprinted in appendix IB. II 
For T > 0, the RHS of the flow equation ()4.30j) is more complicated: a matrix product 
of an inverse tridiagonal matrix (which is a full matrix), a tridiagonal matrix Q, and 
another inverse tridiagonal matrix. In a forthcoming article and in appendix IB.2| I 
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present an algorithm to compute even this product of full N x N matrices in O(N) 
time |Enss 2005j . However, the RHS, which has jumps in A at every Matsubara 
frequency (see section l4.2.3j) . makes the integration at finite temperature scale as 
0(N/T), such that only systems up to N = 10 4 have been computed so far for low 
temperatures T = 10~ 4 . 

The flow equation (|4.H5jl of the density-response vertex used below for the Friedel 
oscillations at T = is of the same computational complexity as the self-energy 
flow at T > 0: the RHS involves the matrix product of an inverse tridiagonal, a 
tridiagonal, and another inverse tridiagonal matrix. This can also be done in O(N) 
time by the same algorithm in appendix IB.2[ and at T = the density profile rij of 
a system of size N = 10 7 can be computed in a few hours for each j. 



4.3 Flow of the density-response vertex 

The expectation value of the local density rij could be computed from the local 
one-particle propagator Gjj if G was known exactly However, the approximate flow 
equations for E can be expected to describe the asymptotic behavior of G correctly 
only at long distances between creation and annihilation operators in time and/or 
space, while in the local density operator time and space variables coincide. In the 
standard RG terminology rij is a composite operator, which has to be renormalized 
separately. 

To derive a flow equation for rij, we follow the usual procedure for the renormal- 
ization of correlation functions involving composite operators: one adds a term (fijUj 
with a small field <fij to the Hamiltonian and takes derivatives with respect to 4>j in 
the flow equations. The local density is given by 

n,A = — 




where £l(4>j) is the grand canonical potential of the system in the presence of the field 
<f)j. Note that I use the same symbol rij for the density operator and its expectation 
value. In the presence of a cutoff A the grand canonical potential obeys the exact 
flow equation 

d A n A = T tr [e itj0+ [d A Q A (iLu)] [G a (iuj) - G a (ilu)]} , 

which follows from the flow equations for the vertex functions, equation (|2.30jl . and 
the relation between the grand canonical potential and the zero-particle vertex, Q A = 
—Jq- At zero temperature (which is the only case I consider here) the Matsubara 
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frequency sum becomes an integral which, for the sharp frequency cutoff (|4.22jl . can 
be carried out analytically. This yields 

d A n A = ^J2 tT { e "° + M 1 ~ sA M] } > ( 4 - 33 ) 

LU = ±A 

which is the first term of equation (|2.34|) . Because at T = any perturbation of 
the Hamiltonian can be shifted between the bare propagator and the self energy, in 
this section I choose to attribute <pj to the interaction part of the Hamiltonian, not 
to Hq, such that Go remains independent of <fij. The self energy is modified via the 
additional local and frequency-independent contribution (pj Sjjr to its initial value 
Ef?, at scale A . 

The density profile can be obtained from the above equations and the flow equation 
for S A by computing the shift of Q A generated by a small finite perturbation (f>j, i.e., 
by numerical differentiation (cf. section l3.3.3j) . Alternatively, one may carry out the 
4>j derivative analytically in the flow equations, which yields a flow equation for the 
density in terms of the density response vertex. Taking the (f)j derivative in equation 
(HH yields 

8 A nf = ~^J2 tT [e iM+ G\iuo) Rf(iuj)] (4.34) 

LO=±A 

with the density-response vertex 

, dZ A (iuj) 
R A iu) = — -A_^ 



and the propagator G A as defined in equation (j4.23|) . i.e., in the absence of 4>j- We 
compute the self energy S A in the presence of (f>j within the same approximation as 
previously. It is thus determined from the flow equation (|4.24|) with a frequency- 
independent two-particle vertex T A . Taking a derivative of that equation with respect 
to <pj at (f)j = yields the flow equation for the response vertex 

ui=±K 2,2' 3,3' 

Note that R A is frequency independent in our approximation and that there is no 
contribution from the (fij derivative of T A since we neglect self-energy contributions 
in the flow of the two-particle vertex. 

For spinless fermions with a (renormalized) nearest-neighbor interaction, the ma- 
trix R A is tridiagonal, i.e., only the components R A -u and R A . ll±1 are non-zero, and 
their flow is given by 

d ^ R f;l,l = ~1^Z ^2 Y2 Y2 Y2 U tl+r G t+r,l'( iu} ) R f;l',l'+r> G ^+r' ,l+r( iu ) 
u;=±A I' r=±l r'=0,±l 
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9aR%, 1±1 = - H E M'M^w^M' ( 4 - 35 ) 

w=±A /' r'=0,±l 

Then, the flow of is 

dAnf = —— E E E ^+r',l'( iuJ ) R f;l',l'+r' ■ 
lu=±A I' r'=0,±l 

Although naively the flow of R A scales as 0(N 2 ) in time because of the unrestricted 
loops over I, I', it can be performed in O(N) by the method described in the appendix 
IB.21 The initial condition for the response vertex is Rf°ui = o~jio~w- The initial 
condition for the density is = |, for any filling, due to the slow convergence of 
the flow equation (|4.H4jl at large frequencies, which yields a finite contribution to 
the integrated flow from A = oo to Ao for arbitrarily large finite Ao, as in the case 
of the self energy S A °. 

To avoid the interference of Friedel oscillations emerging from the impurity or 
one boundary with those coming from the (other) boundaries of the system one 
suppresses the influence of the latter by coupling the finite chain to semi-infinite 
non-interacting leads, with a smooth decay of the interaction at the contacts (cf. 
section 14. 1 . 1ft . 

4.4 Computation of the conductance 

The linear-response conductance is defined via the infinitesimal current induced by 
an infinitesimal voltage drop at zero bias voltage, G = dl/dV. It is a global quan- 
tity defined over the whole wire, from one lead through the interacting (scattering) 
region to the other lead. Even for a perfectly clean wire (interacting or not), the 
conductance is limited to e 2 /h per channel. For a Id system of spinless fermions 
there is exactly one channel, while for spin-| fermions there are two channels. 

I choose the conductance as the appropriate observable for transport, as opposed 
to the local quantity conductivity, because I am interested not in bulk properties 
but the effect of a specific spatial setup of impurities at defined positions, the double 
barrier, for which there are experimental data available |ptygd 2UUlj . 

4.4.1 Kubo formula 

The conductance is computed, just as the linear-response conductivity, via the Kubo 
formula from the current- current correlation function, see for instance |Mahan 2000} 
chapters 3 and 7]. This correlation function at finite frequency and temperature is 
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defined as 

rP 



jx(iuj) := / dre^ T (T T J R (T)J L (0)} 
Jo 



with Jl,r the current operators at the left and right ends of the system. The retarded 
correlation function 7r rct (^) is obtained by analytical continuation iu i— > u + iO, and 
the dc conductance is given by finally taking the limit u — > 0, 

G,= % lim . (4.36) 

It is important that one does not set u = from the beginning because the limits 
uj — > and q — > (macroscopic transport from one end of the system to the other) 
do not commute |Luttinger" T964| . Another way to compute the conductance at 
zero temperature via persistent currents is explained in |MedenfcSchollwock 2003a| 
IMedenfcSchollw5ck 2003lo| . 

Consider the following microscopic setup (cf. figure l4~Tl on page lo"^ . The scattering 
region on the lattice sites 1, . . . , N ("sample") is interacting; it is connected to semi- 
infinite, non-interacting leads at the interfaces at sites 1 and N. The current and 
total number operators are 

J L := it L (c\c - cjcj 

Jr '■= ^i?( c jv+l c Af — c jv C V+l) . A 

(4.37) 



N 



n c 



with Jl the current flowing from the left lead into the sample, Jr the current flowing 
out from the sample into the right lead, and nc the total particle number in the 



interacting region. 3 



In the exact theory, the current- current correlation function can be expressed in 
terms of the bare current operators Jl,r and the Ipi two-particle interaction vertex 
T, see figure l4~3l The expressions for these two contributions to n{iuj) are 



TT^Hiuj) 



-(it L )(it R )T ^ [Gjv,i(ie + iu) G 0>N+1 (ie) + G N+lfi (ie + iu) G ltN (ie) 



3 Thc components of the gauge potential = (<p(t,j),A(t, + 1})) are the scalar potential ip 
at time t and on site j, and the vector potential A on the (directed) bond from site j to site 
j + 1, respectively, where we have introduced the notation [j, j + 1] for the bond. Coupling to 
the potential A^ changes the bare dispersion £ to [Zinn- Ju stin 2002] chapter 34] 



= -A* - e<p(t,j) 

Derivatives with respect to a A([0, 1]), A([N, N + l]) and a global if yield the expressions l|4.H7|l . 
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18+IV 




'R 



Figure 4.3: Two contributions to the current- current correlation function 7r(iu>). The 
shaded region represents the vertex part Ti^z^ie, ie'; ioS). 



G N+lj i(ie + iu) G 0>N (ie) - G N , (ie + iu) G ltN+1 {ie) 



N 



7T (b) (iuj) 



-(^ L )(^)T 2 ^ £ 

ie,ie' Ji,...,j4=l 

G jul (ie + iu) G j 4 (ie) - G jlt0 (ie + iu) G ljh (ie) 
G Nj2 (ie' + iuj) G j3jN+1 (ie') - G N+1>j2 (ie' + iu) G j3>N (ie') 



Kubo formula projected into the interacting region 

Using equations ()4.5j) . f)4.6|) . and ()4.7j) . we can express Ti^ a,b \iu) in terms of Green 
functions with indices in the interacting region 1, . . . , N, 



T Ai(ie, ie + iu) Gi^fa) A/?(ie, ie + iu) Gn,i{^ + iw) (4.38) 



A' 



7T {b) (lUJ 



iu) = T^2 ^2 ^{ie,ie + iu) Gi^{ie) P R;j4tjl (ie,ie + iu) G jltl (ie + iu) 

31,34=1 

with bare current vertices on sites 1 and N, respectively, 

\ L (ie, ie + iu) = -it 2 L [g L (ie + iu) - g L (ie)\ 
\ R (ie, ie + iu) = +it 2 R [g R (ie + iu) - g R (ie)} 

and the current-vertex correction 



N 

ie' 32,33=1 

x Gj 3) Ar(ie / ) X R (ie'; ie' + iu) GNj 2 (ie' + iu) 

and likewise for Pl. 
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Analytical continuation of the Kubo formula 

The Green function ()4.3|) of the lead has a branch cut at the real axis, 

g LtR (e ±i5) = -(e + /i±i5 T i\J^- (e + \i ± i$) 2 ) 
= i( e + /i T V4-(e + /i) 2 )+0(5). 

Thus, for small a;, 

A i!jR (e±ia;,e±iO) = 
A L (e + iu, e - z'0) = +t£ ^4 - (e + + 
A R (e + zw, e - z'0) = -t% - (e + + £>(u/), 

and only the current vertices with incoming and outgoing legs on opposite sides of 
the branch cut contribute in the limit u — > 0. The Matsubara sum in the bubble 
term (|4.38|) is evaluated by a contour integral as explained in |Mahan 20001 chapter 

7], 

x ( + A L (e + z'0, e + iu) Gi >N (e + z'0) \ R (e + z'0, e + iu) G NA (e + iu) 

— A^(e — z'0, e + zcj) Gi^(e — z'0) A#(e — z'0, e + z'u;) Gn,i(c + iw) 
+ (ji,iv(e — ioj) A/j(e — z'u;, e + iO) Gi^e + iO) Ax(e — iu, e + iO) 

— Gi i Ar(e — za;) A#(e — iu, e — iO) G Xi jv(e — iO) A^(e — iu, e — z'0)) . 

The first and fourth term on the RHS have frequency arguments on the same side 
of the branch cut and are, therefore, of 0(u 2 ), hence they will vanish in the limit 
u — ► in equation ()4.3fij) . We retain the other two terms and perform the analytical 
continuation of the external frequency, iu i— > u + iO, 

x ( - A L (e - z'0, e + u + z'0) G 1)N {e - iO)X R (e - z'0, e + u + zO) Civ,i( e + ^ + *0) 
+ A L (e - u - iO, e + z'0) Gi,7v(e - w - ?0)A fi (e - u - iO, e + z'0) Gjv,i(e + iO)) . 

We then substitute e i— > e + in the second term, 

*S(«0=< / ^[/(c + w)-/(c)] 
x A L (e - z'0, e + u + z'0) G 1)JV (e - z'0) A^(e - z'0, e + u + z'0) Gjv,i(e + + zO) . 
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4-4 Computation of the conductance 



As the final step we take the limit uo — > 



G {a) = - lim 




(e, T) is the transmission probability at temperature T without vertex correc- 
tions, 



where we have used G\^{e — iO) GV,i( e + *0) = |GV,i(e + i0)\ 2 : because the Hamil- 
tonian is time- reversal invariant, the amplitude from site 1 to N is the same as from 
N to 1, and G^y is symmetric (not hermitean). The energy integration extends over 
the band of the non-interacting leads, while the factor [4 — (e + /i) 2 ] from the density 
of states of the leads suppresses the transmission towards the edge of the band. 

The f RG provides an approximation of the frequency- independent self energy 
J2jji(T) at zero or finite temperature. S acts as an effective static potential by 
which non-interacting electrons are scattered. The full propagator is determined via 
the Dyson equation (|4.10|) . 

4.4.2 Vertex corrections 

The second contribution (e, T) to the transmission is due to current- vertex cor- 
rections. It is obtained from ir^ following the same steps as for T^ a \ We have 
approximated the full effective two-particle interaction vertex T by a renormalized 
nearest-neighbor interaction with all external frequencies set to zero, i.e., without 
branch cuts. Then the vertex corrections are, omitting the lattice indices and de- 
noting the loop summation by the trace, 



rW(e, T) := t\ t\ [4 - (e + /i) 2 ] \G N}l (e + t0)\ 2 



(4.40) 




- G(e - iO) X R (e - iO, e' + iu) G(e + iuj) Y 
+ G(e' - iu) X R (e' - iuj, e' + iO) G(e + iO) T 

- G(d - iu) X R (e' - iu, e' - iO) G(e' - iO) T) 
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4 Functional RG technique in one dimension 



Again, the first and fourth term on the RHS are by an 0(u) smaller than the other 
two, so we retain only the second and third term and perform the analytical contin- 
uation iu i — > u -\- iO. Substituting e' i— > e' + u in the second term, 



P R (u + iO)=i [ ^[f( e ' + u)-f(e')} 

J 27T^ v ' 



<D(w) 

x tr (G((? - iO) X R {e' - iO, e' + u + iO) G(e' + + z'O) r) ^> 

vanishes since there is no division by u as in equation ()4.36|) . Because T is frequency 
independent, there are no vertex corrections, hence equation (|4.40|) is the complete 
transmission probability in our approximation. 

Conformance with Ward identities 

The approximation that T is frequency independent has another consequence: by the 
flow equation, it follows that ImS = 0, i.e., we do not capture inelastic processes 
at second order in the interaction. They could be included in the flow equation by 
retaining the frequency dependence and imaginary part of V. 

However, the fact that the vertex corrections and ImS vanish simultaneously 
shows that our approximation is at least consistent with the (non-perturbative) Ward 
identity associated with global particle number (charge) conservation. The global 
continuity equation for the interacting region is dnc/dt + J R — Jl = 0. Following 
Oguri 2001 , we define the number and current response functions as the time- 



ordered expectation values of the current and number operators ()4.37|) with two 
extra electron legs (1 < j,j' < N), 

^C;jj'(r;r 1 ,T 2 ) = (T T [n c (r) - (n c )] c^n) ct,(r 2 )) 
$w(t;ti,t 2 ) = (T t J^t) ^(tOcJ,^)) 
*/?;ji'( r ; r i,i" 2 ) = (T t J r (t) cjin) ct,(r 2 )) . 

Performing a Fourier transform on each of these {x = L, R, C), 

®x(t; n, r 2 ) =T 2 J2 ie + iu) e~ ie{T1 ~ T) e -^+")(^). 



IC.lUl 



The Ipi response vertices are obtained by amputating full propagators, 
A x (ie, ie + iu) := [G(ie)] _1 $ x (ie, ie + iu) [G(ie + iu)}' 1 

which can be written in terms of the bare current and the vertex corrections as 

A L n'(ie, ie + iu) = A^fie, ie + iu) 5^ i<L 1 + Pl ivfie, it + ioS) 

>jj ' JJ (4-41) 

A R jj/(ie, ie + iu) = X R (ie, ie + iu) Sf^Sj^N + PR,jj'(ic, ^ + w^) • 
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4-4 Computation of the conductance 



The Ipi Ward identity (jH.12j) in terms of these response vertices then reads 

iuAc(ie, ie + icu) + zA#(ze, ie + icu) — iAi(ie, ie + icu) 
= G _1 (ie + zw) -G~\ie). 

Clearly, the knowledge of the RHS is not sufficient to determine Ac, Ar, and 
separately but fixes only the difference. We use the Dyson equation f)4. on the 
RHS, continue analytically to ie + iu i— > e + uj + iO and ie i— > e — zO, and take the 
limit uj — > 0. Thereby, the density-response term is suppressed, and we are left 
with a relation between current and self energy, using equation ()4.8|) for the lead 
contribution, 

A R (e - z'0, e + iO) - A L (e - iO, e + z'0) 

= -Im£(e + z0) - Im S lcads (e + iO) +ImS(e-iO) + Im E lcads (e - iO) 
= X R \N)(N\ - A L |1)(1| - 2ImS(e + zO) . 

By equation (IPT1) . 

P^(e - iO, e + zO) -P L (e- zO, e + zO) = -2 Im E(e + zO) . (4.42) 

Thus, the current-vertex corrections are related to the imaginary part of the self 
energy. At T = 0, both sides vanish exactly because there is no inelastic scattering, 
while at higher temperatures, neglecting inelastic processes is a consistent approx- 
imation we made. If one tried to improve the approximation by computing vertex 
corrections but keeping a real E, it might be no improvement at all because one 
would violate number conservation, as explained in |BaymfcKadanoff 1961| . 

4.4.3 Different flow schemes 

In order to check the robustness of our fRG results, I have considered also the 
temperature-flow scheme (cf. section 14.2.4)1 where the temperature is successively 
lowered during the fRG flow. Another option is the recently introduced interaction- 
flow scheme (cf. section r4.2.5|) which slowly switches on the interaction strength dur- 
ing the flow (at finite temperature, or even at zero temperature because in the model 
at hand, the finite size already provides a regularization). Using these schemes, I 
have obtained the same results as for the Matsubara-frequency cutoff presented in 
section H.2.HI for a few test cases (cf. figure l4~4~j) . The frequency cutoff scheme at 
finite temperature remains numerically the most efficient method. 
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10" 4 10" 3 10" 2 10" 1 10° 10 1 10 2 
T 

Figure 4.4: Comparison of frequency cutoff and temperature flow /RG schemes to first- 
order perturbation theory for the resonant conductance through a double bar- 
rier (N = 10 4 , U = 0.5, V = 10, N dot = WO). This shows that different 
/RG schemes essentially give the same results, indicating the robustness of the 
method. 
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5 Luttinger liquids with impurities 



The behavior of three-dimensional metals can be described by Fermi-liquid the- 
ory, where the elementary excitations are fermionic quasi-particles. In one di- 
mension, however, the situation is completely different (for reviews see |Voit 1995j 
IGiamarchi 20fT3j). 

A model for interacting Id spinless fermions with linear dispersion was introduced 
by |Tomonaga 1950| and |Luttinger 1963 and has been solved by [MattisfcLieh 1965j . 



It is characterized by a gapless spectrum of collective bosonic, density-wave elemen- 
tary excitations with a linear dispersion, as well as charge and current excitations, 
and correlation functions decaying at zero temperature algebraically in space and 
time. |LutherfcPeschel 1974aj and |Mattis 1974] introduced the technique of boso- 
nization as a particularly simple method to describe the low-energy properties of 
this model. Any non-linearity of the fermion dispersion, for instance on the lattice, 
gives rise to interactions between the collective modes. For spinless fermions, a per- 
turbative expansion of these bosonic couplings is completely regular in the infrared, 
and the low-energy excitation spectrum is similar to the Luttinger-model spectrum 
|Haldane 1981] . Therefore, the term Tomonag a- Luttinger liquid (tll) was coined for 
the generic low-energy phase of interacting Id fermions. The correlation functions 
obey power laws with exponents that are functions of a single interaction parameter 
K with K = 1 in the non-interacting case, < K < 1 for repulsive and K > 1 for 
attractive interaction. 

In the presence of impurities or boundaries there are further characteristic power 
laws. For example, the local density of states Pj(cu) = — ^lmGjj(uj + iO) near a single 
impurity or boundary is suppressed for repulsive interaction as [KanefcFisher 1992cj 

pj{uj) ~ \u)\ aB for lu — > 

where 

a B = l/K - 1 

is the boundary exponent. The linear conductance across an impurity of arbitrary 
strength scales as [KanefcFisher 1992aj 

G(T) ~ T 2aB for T -> 0. 
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5 Luttinger liquids with impurities 



For repulsive interaction (ocb > 0), a system containing even a very weak impurity 
is insulating at T = 0, effectively cutting the chain into two disconnected parts. 
Another TLL signature is the decay of Friedel oscillations: an impurity or boundary 
induces oscillations in the density profile rij whose amplitude scales with the distance 
EggerfcGrabert 1995| 



x as 



Arij ~ x K for x — ► oo, 

in contrast to Arij ~ x~ l for a Id Fermi gas. These power laws are strictly valid 
only in the low-energy limit. This raises the question at which scales N, T the 
asymptotic TLL behavior sets in for a specific model, and how the system behaves 
before reaching the asymptotic limit. 

One way to study TLL behavior is to consider the Id model of spinless fermions 
on the lattice with nearest-neighbor interaction U, which we have defined already in 
section E~T1 Without impurity, this model can be solved exactly by the Bethe ansatz 
Yangfc Yang 1966] ; it is a Luttinger liquid at any interaction strength U and any 



filling, except for \U\ > 2 at half filling. The TLL interaction parameter K is given 
for \U\ < 2 at half filling by jHaldane 1980] 

K = — arccos 

7T V 2 

For U > 2 at half filling, a phase transition towards a charge- density wave occurs, 
while for U < —2 the system undergoes phase separation. In this work we shall 
concentrate on TLL physics, so parameters have to be chosen to stay away from 
these phase transitions. 

The local density of states has been studied in depth by |Andergassen 2005 ; I 



will present new results for Friedel oscillations (section 15.1)1 and transport through 
double barriers (section 15.2)1 . The fRG turns out to be a versatile tool to study 
physical effects on energy scales ranging over several orders of magnitude, at weak 
to intermediate interaction strength. 



5.1 Friedel oscillations 

In a normal Fermi-liquid metal, impurities induce Friedel oscillations in the density 
profile jFriedel 1958) ITuttofeZawadowski 1985) which far away from the impurity 
have the form 

An(x) ~ cos(2A; F x + 5) x~ d 

where An(x) = (n(x) — no) /no is the normalized density profile, fcp is the Fermi 
wave vector, x the distance from the impurity, d the dimension of space, and 5 a 
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5.1 Friedel oscillations 



phase shift. This raises the question which density profile is generated in a Luttinger 
liquid. 

The continuum Id Luttinger model with a single impurity is integrable, hence 
Friedel oscillations should be computable exactly for any coupling strength. However, 
this is technically difficult, and only approximate results have been obtained except 
for an exact solution at K — ^ |lls 19 96 . The asymptotic behavior for very weak 
and very strong impurities was studied by jKanefeFisher 1992al IKanefcFisher 1992bj 
using the RG method described below. On the other hand, the crossover at interme- 
diate impurity strength is particularly important for understanding transport in Id 
wires and was analyzed for weak interaction by |myo 199*3] . 

|EggerfcGrabert 1995| first studied arbitrary impurities for repulsive interaction 
using bosonization, which is valid at low temperatures. For a strong scatterer, the 
amplitude of the Friedel oscillations decays as x~ K . For a weak scatterer, there is 
a crossover from the asymptotic x~ K decay for large distances (x ^> xq) to linear- 
response decay |Voit 199*5] as x l ~ 2K for short distances (x -C xq)- The crossover 
scale diverges as x ~ V^ 1 ^ 1 '^ for V — > 0. 

Shortly after, |lls 1995] obtained exact results for the density profile in the con- 
tinuum for K — | at arbitrary temperature. For this particular interaction, the 
problem can be mapped to a free-fermion model, which simplifies the calculation. 
At T = 0, 

An{x) = Ahcos(2k F x + r? F ) e** h '*K {torh 2 x) 

with impurity strength h = X/y/2 (A is the coefficient of the cos[0(O)] term in the 
bosonized Hamiltonian) , K (x) is a modified Bessel function and rjp = —KirX/kp 
a phase shift. For large distances, the amplitude decays like x~? as expected. At 
T > 0, 

A . , , , , / 4vrT (\ 1 h 2 1 - coth(27rTx)\ 

An(x) = cos(2/c F x + Vf) \ — r~, — ^r^ F -, -; 1 + 2—, 

v ' v ' ; y sinh(27rTx) V 2 2 T 2 / 

with the hypergeometric function F(a,b;c,d). For x ^> 1/T the amplitude of the 
Friedel oscillations decays exponentially. 

Recently, |gyl 2UU"ij used the functional bosonization technique and a self-consistent 
harmonic approximation at low temperature and weak impurity strength to obtain 
the Green function at arbitrary interaction K. It shows the full crossover from the 
impurity-dominated behavior at short distances to the pure tll behavior at large 
distances. From the Green function, one can extract the density profile as 

An(x) = cos(2A;fx + t]f) \smh(2nTxK/vF)\~ K . 



81 



5 Luttinger liquids with impurities 



For distances shorter than the thermal coherence length, x <C ttvf/T, the amplitude 
decreases as x~ K , while for larger distances it is suppressed exponentially. 

In section EOl above I have described how to compute the density profile in the fRG 
framework, illustrating the renormalization of a composite operator and the flow of 
response functions in a concrete example. In the following section I will present my 
results for the density profile for a wide range of parameters and compare them with 
exact dmrg results and the asymptotic formulae from bosonization presented here. 

5.1.1 Results 

Figure 15.11 shows fRG and dmrg results for the density profile rij for a spinless- 
fermion chain with 128 sites and interaction strength U = 1 at half filling. The 
Friedel oscillations emerge from both boundaries and interfere in the center of the 
chain. The accuracy of the fRG results is excellent for all j. 

For incommensurate filling factors the density profile looks more complicated. This 
can be seen in figure 15.21 where fRG results are shown for the density modulation 
\rij — n\ near the boundary of a system with an average density n = 0.393 and 8192 
sites. For long distances from the boundary the oscillation amplitude has a well- 
defined envelope which fits to a power law as a function of j. In the following I will 
examine the large-distance behavior of the amplitudes more closely for the half-filled 
case. 

Figure 15.31 shows fRG results for the amplitude of density oscillations emerging 
from an open boundary, for a very long spinless fermion chain with 2 19 + 1 sites and 
various interaction strengths U at half filling. The other end of the chain (opposite 
to the open boundary) is smoothly connected to a non-interacting lead. In a log- 
log plot (upper panel of figure 15. 3 j) the amplitude follows a straight line for almost 
all j, corresponding to a power-law dependence. Deviations from a perfect power 
law can be seen more clearly by plotting the effective exponent atj, defined as the 
negative logarithmic derivative of the amplitude with respect to j (see the lower 
panel of figure l5~3|) . The effective exponent is almost constant except at very short 
distances or when j approaches the opposite end of the interacting chain, which is 
not surprising. From a comparison with the exact exponent (horizontal lines in the 
figure) one can assess the quantitative accuracy of the fRG results. 

Effective exponents describing the decay of Friedel oscillations generated by site 
impurities of various strengths are shown in figure 15.41 for a half-filled spinless 
fermion chain with 2 18 + 1 sites and interaction U — 1. Both ends of the inter- 
acting chain are coupled to non-interacting leads to suppress oscillations otherwise 
induced by the boundaries. For strong impurities the results are close to the bound- 
ary result (cf. figure l5~^|l . as expected. For weaker impurities the oscillations decay 
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more slowly, i.e., with a smaller exponent, and do not reach the boundary behav- 
ior within the range of our chain for V < 1. For very weak impurities (V = 0.01 
in figure I5.4J1 the oscillation amplitude follows a power law corresponding to the 
linear-response behavior with exponent 2K — 1 at intermediate distances. 

The same crossover between linear-response behavior for very weak impurities and 
tll behavior for strong renormalized impurities is observed in the oscillations of the 
effective impurity potential E^- (cf. figure I5~K|) . When the density oscillations decay 
with exponent K, the oscillations decay with exponent 1. In the linear-response 
regime, however, the respective exponents are 2K — 1 and K. This was explained in 
jMMSS 2002al IMMSS 2002bj: 




Figure 5.1: Density profile nj for a spinless fermion chain with 128 sites and interaction 
strength U = 1 at half filling, /rg results show an excellent agreement with 
numerically exact dmrg data. 



since the linear-response backscattering amplitude scales with the cutoff A as Vk Y -k ¥ ~ 
(1/A) 1 ~ K , this provides the generic scaling law in the low-energy regime, and also the 
self energy scales away from 2/cp with the same exponent, S^^' ~ {k — k' — 2kpY^ K . 
Performing a Fourier transform, the real-space decay scales as \j — jo\~ K ■ 

Finally, I present results for the effective exponent of the density-oscillation decay 
in the case of an attractive interaction U — —1, see figure 15.61 In that case the 
effective impurity strength should scale to zero at low energies and long distances 
|KanefcFisher 1992aj . Indeed, for weak and moderate bare impurity potentials the 
effective exponent in figure IB~6l approaches the linear-response exponent 2/^—1. Only 
for very strong impurities the density oscillations decay with the smaller exponent 
K over several orders of magnitude. 
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Figure 5.2: Density modulation \rij — n\ as a function of the distance from a boundary, 
for spinless fermions with interaction strength U = 1 and average density 
n = 0.393 on a chain with 8192 sites; the dashed line is a power-law fit to the 
envelope of the oscillation amplitudes with exponent K = 0.785 corresponding 
to as ~ 0.274. It shows that Friedel oscillations display the expected behavior 
also away from half filling. 



5.2 Transport through double barriers 

The motivation to study double barriers using the fRG is twofold: {%) the double 
barrier is physically interesting and experimentally accessible, for example in specifi- 
cally fabricated quantum wires or double kinks in carbon nanotubes, see for instance 
|ChamonfcWen 1993[Iptygd 200T] . and (ii) it is a multi-scale problem; in addition 
to the temperature T, there are additional scales such as the size of the dot region 
Ndot between the barriers and the strength of the barriers Vl,r, cf. figure lo~71 The 
double barrier exhibits universal scaling in certain limits — we observe several differ- 
ent power laws — but also non-universal crossover behavior in between, as different 
physical processes become relevant. This is particularly important because in exper- 
imental setups, the parameters are usually in the intermediate range. Therefore, the 
double barrier provides a showcase for the power of the fRG to treat all scales on an 
equal footing. Nevertheless, as a note of caution, several important ingredients for 
a realistic description are still missing in our model, although they can in principle 
be handled by our method, such as the spin degree of freedom, higher-dimensional 
leads, and realistic contacts. We therefore refrain from a detailed comparison of our 
findings to experiments. 

For a weak single impurity, the conductance in the limit T — > is suppressed for 
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10 1 10 2 10 3 10 4 10 5 

j 



Figure 5.3: Amplitude (envelope) of oscillations of the density profile rij induced by 
a boundary as a function of the distance from the boundary, for spinless 
fermions with various interaction strengths U at half filling; the interacting 
chain with 2 19 + 1 sites is coupled to a semi-infinite non-interacting lead at 
the end opposite to the boundary. Upper panel: log-log plot of the amplitude. 
Lower panel: effective exponents for the decay, and the exact asymptotic ex- 
ponents (Bethe ansatz) as horizontal lines. The plateaux allow to read off the 
exponents very accurately. 



repulsive interaction and enhanced for attractive interaction as compared to the non- 
interactine; case jLutherfcPescfiel 1 974bl|ApelfeRice 19821 - jKanefcFisher 1 992aj showed 
that this is a universal result independent of the impurity strength: for repulsive 
interaction, the conductance across the impurity scales as G(T) ~ T 2 ° B , such that 
asymptotically for T — > the wire becomes insulating. This was derived using an RG 
method: a weak barrier, or rather its backscattering component V(2k-p), is a relevant 
perturbation of the clean system and grows stronger in the RG flow. On the other 
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Figure 5.4: Effective exponent for the decay of density oscillations as a function of the 
distance from a site impurity of strengths V = 0.01, 0.1, 0.3, 1, 10 (from 
bottom to top); the impurity is situated at the center of a spinless fermion 
chain with 2 18 + 1 sites and interaction strength U = 1 at half filling; the 
interacting chain is coupled to semi-infinite non-interacting leads at both ends. 
This shows the crossover from the density response regime (V <C 1) to the 
boundary behavior (V 3> 1). 

hand, a weak link between two otherwise separate semi-infinite chains is an irrele- 
vant perturbation which remains weak in the RG flow. Because in the weak-impurity 
and weak-link limits the direction of the RG flow is compatible towards a strong 
barrier or weak link, respectively, Kane and Fisher connected both perturbatively 
accessible limits and concluded that any barrier becomes strong in the asymptotic 
limit. This is supported by an exact solution at K — |. Conversely, an attractive 
interaction suppresses an initial backscattering to yield perfect transmission even for 
a system with an impurity. The fRG has been used successfully to reproduce the 
one-parameter scaling of the conductance for 1/2 < K < 1 |mamss 2 003 . 

In the case of a double barrier, there are resonance peaks in the conductance 
G(Vc) as the gate voltage Vq is varied. Using an RG analysis similar to the one for 
the single barrier for asymptotically low temperatures, [KanefcFisher 1992cj derived 
a phase diagram for the peak conductance G p depending on the barrier strength 
V and interaction parameter K. For attractive interaction K > 1 there is perfect 
transmission for any Vq, while for repulsive | < K < 1 there are sharp conductance 
peaks at which perfect transmission is reached. For K < | no transmission is possible 
for any Vq, and for | < K < | there is a line (Kosterlitz-Thouless separatrix) of 
critical barrier strengths: for stronger barriers no transmission is possible, while for 
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Figure 5.5: Effective exponent for the decay of oscillations of Hj t j as a function of the 
distance from a site impurity of strengths V = 0.01, 0.1, 0.3, 1, 10 (from 
bottom to top), for the spinless fermion model at half filling and interaction 
strength [7 = 1; the impurity is situated at the center of a chain with L = 2 18 +1 
sites. 



£= 
CD 
£= 
O 
CL 
X 
CD 



1.8 



1.6 



1.4 



1.2 




-i-iii'ffS'SS'S'S □ □ □ ^ 

El. 



10 1 10 2 10 3 10 4 10 5 



J-Jo 



Figure 5.6: Effective exponent for the decay of density oscillations as a function of the 
distance from a site impurity of strengths V = 0.1, 1, 10, 100, 1000 (from 
top to bottom) for the same chain as in figure \5J\ but now with an attractive 
interaction, U = — 1. This demonstrates that our method works well also for 
negative U . 



weaker barriers conductance peaks appear. The line shape of the resonance peaks is 
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Figure 5.7: Model of the double barrier. 



non-Lorentzian, the tails falling off as (AVq) 2 ^ k , where AVq is the detuning from 



resonance. 



Shortly afterwards, |FurusakifcNagaosa 1993| studied the double barrier for arbi- 
trary interaction by second-order perturbation theory in the barrier strength V <C 1. 
At low temperatures, the deviation of the resonant G P (T) from e 2 /h scales as T 2K . 
As the temperature is increased, perturbation theory in the inverse barrier strength 
1/V <C 1 and a Master-equation approach [Furusaki 1998j yield a regime of uncor- 
rected sequential tunneling (ust) characterized by peaks of height G P (T) ~ j^b-i 
and width w(T) ~ T, for T <C A dot , where A dot := nv F /N dot is the level spacing of 
the dot. For Adot <C T <C B, where B is the bandwidth, G P (T) increases as T 2aB 
for increasing T. 

However, recent experiments on carbon nanotubes |ptygd 2UUlj have reported 
for the temperature range r < T « Adot that both the resonant peak height 
G P (T) ~ T - 7 and G(Vg) peak width w(T) ~ T decrease with decreasing tempera- 
ture. Here, V is the width of a resonance in the transmission probability T(e) given 
by equation (|4.40|) . which for strong barriers is related to the tunneling rate into and 
out of the dot. While also |FurusakifcNagaosa 1993 iFurusaki 1998] found power laws 
with a positive power of T, they were expected to occur only for very strong repulsive 
interactions, where olb is at least twice as large as the value as ~ 0.6 . . . 1.0 observed 
in carbon nanotubes. |ptygd 2UUlj suggested that correlated sequential tunneling 
(cst) with exponent — 1 > dominates over UST with exponent as — 1 < 
in the temperature range r < T < A dot . Subsequently, jTGCPD 2UU2l Iteg 2UU4j 
argued that cotunneling processes of second order in the end-tunneling local density 
of states (each with exponent as) are the leading contribution to the peak conduc- 
tance (hence the 2a b), and claimed perfect agreement with the experimental data 
of |ptygd 2UUlj . In contrast, [Furusaki 1998j had found that cotunneling via vir- 
tual intermediate states dominated the tails of the conductance peak away from the 
resonance. This led to a renewed interest in transport through double barriers, with 
several new methods available to advance beyond asymptotically low temperatures. 
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Shortly after, |PolyakovfeGornyi 2003 INazarovfcGlazman 2003j studied the dou- 
ble barrier at arbitrary temperature and impurity strength but weak interaction 
using the leading- log resummation of |myg 199*3*1 Iygm 199*""" . They considered the 
RG flow of the energy-dependent scattering amplitudes, based on a resummation of 
perturbation-theory diagrams with leading logarithmic divergences a/ ln n ( | fc _j, F | d ), 
where d is the spatial range of the interaction. Inelastic processes are sub-leading at 
weak interaction. Their results for a single resonant level agreed with |FurusakifcNagaosa 1993| 
for T < A dot , finding UST but not CST. |PolyakovfeGornyi 2003] extended the RG 
study to a multi-level dot and weak barriers (closer to experimental parameters) but 
again confirmed UST. 

Recently, |KomnikfcGogolin 2003| found an exact solution for a special model at 
interaction strength K = |. However, in that model the conductance scales as 
T" 1 even for T < A^t where one would expect UST, as if the fermions were non- 
interacting. Because it shows no sign of either UST or CST, this result is probably 
not generic and cannot help to resolve the puzzle. 

|HiiglefcEgger 2004] computed the conductance using the Quantum Monte Carlo 
(qmc) method for interaction K = 0.6 and obtained data interpreted to be consistent 
with CST for weak barriers and sufficiently high temperatures, which was explained 
by an additional transport channel for strong interaction. 

Against this background, it is desirable to have an unbiased method to compute the 
conductance without restricting from the beginning — by physical intuition — which 
physical processes are dominant. The fRG provides this with the only restriction that 
it is perturbative in the renormalized interaction but contains contributions of all 
orders in the bare interaction U. It confirms the UST picture but in addition allows 
to vary all parameters to see exactly at which temperature, barrier strength, dot 
size etc. the universal scaling sets in, if at all. While the justification of the fRG was 
weak-coupling, we reproduce one-parameter scaling for the single barrier at U — 2 
very accurately |mamss 2003j . and our results for the double barrier qualitatively 
agree with the QMC data (U = y/3), indicating that the fRG is reliable up to this 
interaction strength. 



5.2.1 Results 

The fRG has been used before to compute the conductance through a single im- 
purity at zero temperature |MedenfcSchollwock 2003al IMedenfcSchollwock 2003bl 
Imamss 2003j . After having developed and implemented the finite-temperature fRG, 
I have obtained results for transport through a double barrier, symmetric or asym- 
metric, at or off resonance, with weak or strong barriers, enclosing small or large 
dots, for repulsive interaction (see also |meams 2()()5[ Iemabms 2005] ). I will report 
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my findings for several interesting regions in this large parameter space, showing 
agreement with known results as well as providing clarification of a contentious is- 
sue. 

I first present results for a strong symmetric double barrier, Vl,r = 10, enclosing 
a dot of iVdot = 6 sites. The interacting system size is taken to be N = 10 4 , in 
agreement with the size of realistic samples of carbon nanotubes. As the gate voltage 
Vq is varied across the band, there are six resonance peaks with level spacing roughly 
A dot ( figure 15"%)) . In non- interacting systems the peaks have a finite width even at 
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Figure 5.8: The conductance G(Vg) as a function of gate voltage Vq for N dot = 6, U = 0.5, 
Vl,r = 10, N = 10 4 and different T. One observes N^ot resonance peaks which 
are widened with increasing temperature. 



T = 0, depending on the shape of the barrier and the local density of states of 
the leads at the barrier. Repulsive interaction has a strong influence on the line 
shape of the resonance peaks: they become infinitely sharp in the asymptotic limit 
(N — > oo, T — > 0) but perfect transmission e 2 /h is still possible for symmetric 
barriers. As the temperature is increased the peaks become wider and lower. The 
peak conductance G P (T) as a function of temperature shows several different power 
laws for appropriately chosen dot parameters. All the peaks in figure 15.81 have a 
different shape and size due to band effects, however G P (T) behaves similarly for 
each peak, and in the following I will always consider the peak closest to Vq — 0. 

Figure 1531 shows the peak conductance G P {T) for three different dot sizes N dot in 
the upper panel, while in the lower panel the logarithmic derivative is plotted. For 
temperatures larger than the bandwidth, the conductance scales as G P (T) ~ T _1 . 
This is because /'(e) and the transmission in equation (|4.39j) vary only very little over 
the band, but /'(e) decreases with increasing temperature as \jT . For temperatures 
below the dot level spacing Adot (indicated by the arrows), the peak conductance 
decreases as G P (T) ~ T as_1 (ust). This exponent is marked in the lower panel by 
the dashed line. In the UST regime the conductance integral is dominated by a single 
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Figure 5.9: The resonant peak conductance G P {T) for different dot sizes. Upper panel: 
Ndot = 2 (circles), 6 (squares), and 100 (diamonds), with the respective level 
spacing A dot indicated by the arrows. U = 0.5, N = 10 4 , Vl,r = 10. The solid 
curve shows G{T)/2 for a single barrier. Lower panel: Logarithmic derivative 
of G P (T). Solid line: 2a#/ dashed line: olb — 1; dash-dotted line: las — 1. 
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Figure 5.10: Transmission T(e,T) for several temperatures: left panel T = 0.004, right 
panel T = 0.063. The dashed line is the derivative of the Fermi function 
at the respective temperature. The dot parameters are N^t = 100, U = 1, 
N = 10 4 and Vl,r = 10. This shows that for low temperatures ( USt), a single 
transmission peak determines the conductance, while for higher temperatures 
(Kirchhoff) many peaks contribute. 



peak in the transmission T(e, T) with a width V <T (cf. figure IB~TU1 left panel). As 
the temperature is lowered, the UST regime extends down to T* where T — Y. In 
G(Vg) the peaks are well separated and have a width w(T) ~ T. 

For the above parameters and dots with A^ot > 4, the conductance shows a 
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third regime for Adot ^ T < 1. In this region, the conductance increases with 
increasing temperature, leading to a non-monotonic overall T dependence. For large 
dots iV dot > 30, the conductance increases according to the power law G P (T) ~ T 2oiB , 
indicated by the solid line in the lower panel of figure 15.91 Many peaks in the 
transmission T(e, T) contribute to the conductance (cf. figure l5~TTH right panel). This 
regime is most clearly seen for strong barriers Vl,r 3> 1 and large dots iVdot 3> 1. 

If we add the resistances of both barriers separately according to Kirchhoff 's law, 
1/G = 1/Gl + 1/Gr, we obtain the solid curve in the upper panel of figure l5~Hl As 
our data lie on this curve for T > Adot, we suggest the following interpretation: the 
electrons first tunnel through one barrier with a certain probability, then incoherently 
through the next. Only as the temperature is lowered to T « Adot, the electron 
"sees" both barriers coherently and tunnels resonantly, such that the conductance 
increases towards e 2 /h as T — > 0. 
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Figure 5.11: On- and off-resonance conductance. Upper panel: G(T) for U = 0.5, N = 
10 4 , Vl,r = 10, and N^ot = 100. On resonance AVg = (circles), near 
resonance |AVg| = 0.001 (squares), and in a conductance minimum with 
|AVg| = 0.04 (diamonds). The off-resonance curves lie on the resonance 
curve for T > |AVg|, then cross over to decay as G P (T) ~ T 2 " B . Lower 
panel: Logarithmic derivative ofG(T). Solid line: las', dashed line: as — 1. 
At low temperature the double barrier off resonance behaves like a single 
impurity. 



It is thus only in the low-temperature range where it makes a difference whether 
Vg is tuned to resonance or not (cf. figure l5~TTj) . If the gate voltage is slightly off 
resonance, AVg ^ and \SVg\ Adot, the conductance is identical to the resonant 
case for T > AVg, while for lower temperatures the double barrier acts as a single 
impurity of strength |AVg|, hence the conductance scales as G(T) ~ T 2oiB . 

Up to this point, we have presented conductance data for the weak interaction 
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Figure 5.12: Conductance for different interaction strength. Upper panel: Peak conduc- 
tance G P (T) for U = 0.5 (full symbols) and U = 1 (open symbols). Curves 
are shown for dot sizes N^ot = 2 (red) and N^ot = 100 (blue). Lower panel: 
Logarithmic derivative ofG{T). The dependence of the exponents on the in- 
teraction strength support the claim that the exponents are ag — 1 and 2ob, 
respectively. 



strength U = 0.5 («g = 0.165). The numerical exponent of —0.835 in the UST regime 
was identified with og — 1, and the exponent 0.330 in the Kirchhoff regime with 2a b- 
Figure EH21 shows the exponents for a different value of U: at U — 1 {as = 0.35), 
the exponents —0.65 and 0.70 can be read off, confirming this identification within 
the numerical accuracy. A comparison of the exponents for strong and weak barriers 
can be found in jEMABMS 2005j . 

For strong barriers, large dot size and weak interaction, there is no indication of a 
power law with exponent 2a B -l claimed by |tocpd 2(}Q2|Iteg 207!H|HuglefeEgger 2004| . 
Consider, therefore, a dot with parameters as close as possible to those in |Hiig lcfcEg ger 2004 



weak to intermediate barriers, intermediate dot size and larger U (cf. figure l5~T^j) . 
For a small dot iVdot = 10 (circles), neither the exponent as — 1 (dashed line) nor 2a b 
(solid line) is clearly developed. For larger dots iVdot = 30 (squares) and A^ot = 50 
(diamonds), the as — 1 exponent (ust) is clearly visible, but the Kirchhoff exponent 
2as is still not reached. In the QMC data |HiiglefcEgger 2004| the conductance is 
fitted to a slope of 1/3 in a log-log plot, which is interpreted as scaling with expo- 
nent 2as — 1 (cst) for as = 2/3. Our data shows the conductance curve with a 
slope of 0.3 . . . 0.7, depending on the dot parameters, but no fixed scaling exponent. 
Therefore, we interpret this as non-universal behavior. 

Figure 15. 1 41 shows the resonant peak conductance for strongly asymmetric barriers. 
The Kirchhoff regime is unchanged qualitatively as both barriers are very strong. For 
lower temperatures one observes a UST regime which crosses over into single-impurity 
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Figure 5.13: Weak barriers, Upper panel: G P (T) for U = 1.5, N = 10 4 , Ndot = 10 with 
Vl,r = 0.8 (circles), N^t = 30 with Vl,r = 1-5 (squares), and N dot = 50 mt/i 
Vl,r = 0.8 (diamonds). The arrows indicate A^ oi /or i/ie different dot sizes. 
Lower panel: Logarithmic derivative of G P (T). Solid line: 2as; dashed line: 
as — 1. For small dots, the 2ob power law is not clearly developed; however 
the conductance agrees with the QMC data [Hugle&Eg ger 200$ . 
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Figure 5.14: Strongly asymmetric barriers, Upper panel: G P (T) for Vl = 5, Vr = 50, 
U = 0.5, N = 10 4 , Ndot = 2 (circles) and Ndot = 20 (squares). Lower panel: 
Logarithmic derivative of G p (T). Solid line: 2czb; dashed line: as — 1. 



scaling G P (T) ~ T 2otB for small dots at very low temperatures. 

For a small dot (A^ot = 1) with weak barriers the conductance approaches the per- 
fect value e 2 /h according to the power law [KanefcFisher 1992bl|FurusakifcNagaosa 1993 
(cf. figure l5~To} 
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Figure 5.15: Effective exponent of I — G P (T) / (e 2 / h) for a small dot with N dot = 1, U = 1, 
N = 10 4 , and intermediate to weak hopping barriers ti r . Dash-dotted line: 
2K. 



Thus, our method yields up to four different power laws within the same framework. 



5.3 Summary 

I have demonstrated that the fRG is a powerful method for multi-scale problems. 
As applications I have studied Friedel oscillations off an impurity or boundary in 
an interacting wire, and transport through a double barrier. Both problems were 
previously investigated using effective field-theoretical models; in the cases where 
exact results for such models are known I have obtained quantitative agreement for 
weak to intermediate interactions 1/2 < K < 1. For the Friedel oscillations at 
T = 0, I observed the exponent K at large distances and 2K — 1 for weak impurities 
at short distances jAEMMSS 2UU4j . 

For resonant tunneling, depending on the parameters of the quantum dot, I found 
several temperature regimes with power-law scaling as well as non-universal behavior 
(cf. figure EHEI) |meams 2()()5[ Iemabm s 2005J. All these temperature regimes are 
obtained within the same approximation scheme. The crossover between the regimes 
can be studied in detail. For parameters for which a comparison is possible the results 
agree with the ones obtained in lowest-order perturbation theory in the barrier height 
and inverse barrier height. I did not find any indications of a CST regime with the 
exponent Icxb — 1 predicted from an approximate Master-equation approach, and 
seemingly supported by QMC data. If it were present, our method should be able to 
reveal such a regime since its scaling exponent differs from the UST exponent already 
at leading order in the interaction. 
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Figure 5.16: Schematic plot of the different regimes for the scaling of the conductance 
G P (T) at resonance for symmetric barriers found using our method. The 
upper panels have higher barriers than the lower ones; the right panels have 
larger dot sizes than the left ones. Aw = ttvf/N is the energy scale of the 
interacting wire. 
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In this thesis I investigate the role of symmetries and conservation laws in the func- 
tional renormalization-group formalism, and study specifically Friedel oscillations 
and transport in Id correlated electron systems. 

In the functional renormalization group (fRG) method an energy cutoff scale is 
introduced in the bare propagator as a flow parameter (Chapter |2J. By solving 
coupled differential flow equations for the Green functions, the effective behavior 
on all energy scales can be computed for a given microscopic model. Because the 
full flow-equation hierarchy can be solved exactly only in a few special cases, for 
instance the Luttinger model |sbk 2004] . in most practical applications one has to 
truncate the hierarchy by setting the flow of, say, the three-particle and higher Green 
functions to zero. This approximation is justified perturbatively in the renormalized 
interaction or some other small parameter. For applications such as the 2d Hubbard 
model and Id impurity problems, different basis sets of Green functions have proven 
to be particularly useful. 

The presence of continuous symmetries in the bare action leads by the Noether 
theorem to conservation laws and Ward identities relating Green and response func- 
tions (Chapter El). The solution of the infinite flow-equation hierarchy preserves the 
symmetry once the cutoff is removed. The truncated flow equations with a momen- 
tum cutoff, however, generally violate the Ward identities. After discussing previous 
results from the high-energy physics literature, I show that if a manifestly gauge- 
invariant construction is possible as, for instance, in the temperature-flow scheme, 
the Ward identities between Green and response functions can be satisfied exactly 
despite truncations. 

The related property of self-consistency is satisfied by construction in the con- 
serving approximations |BaymfcKadanofl 1961 , and it would be desirable if it were 
also satisfied in truncated fRG flows. However, I show that the commonly used trun- 
cations generally violate self-consistency. For special reduced models it has been 
shown that truncated flow equations can be modified to satisfy self-consistency and 
even yield the exact mean-field solution |Katanin 20041 Ishml 2004] , so the hope re- 
mains that this may be generalized to arbitrary interaction. However, in the Id 
lattice model of the Luttinger liquid, the truncated fRG is surprisingly successful 
and self-consistency does not appear to play an important role on the level of our 
approximation. 
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Correlated electron systems in one dimension show Luttinger-liquid behavior with a 
strong interplay between interaction and impurities. The fRG method is particularly 
apt to compute the single-particle and transport properties of a lattice model of 
spinless fermions with one or two impurities at all energy scales (Chapter EJ). For 
systems with up to 1000 sites the results agree well with the numerically exact 
dmrg. For large systems the fRG reproduces the asymptotic, universal power laws 
known from thermodynamic Bethe ansatz and bosonization, as well as various field- 
theoretical methods, and in addition indicates the onset of and crossover between 
different power-law regimes. 

I have concentrated on two observables: the spatial density profile for Friedel 
oscillations and the conductance. In the exact solution the density profile can be 
obtained from the full propagator; within our approximation it is much more accu- 
rate to instead treat the density-response vertex as a separately flowing composite 
operator. On the other hand, the conductance is completely determined by the full 
propagator because current-vertex corrections do not arise in our approximation in 
accordance with Ward identities. The practical usefulness of the method relies on 
the fast computation of loop integrals (appendix EI); a little- known mathematical 
trick is used to develop a new algorithm linear instead of quadratic in the system 
size which allows to treat lattices as large as 10 7 sites very accurately. 

The Friedel oscillations off an impurity or boundary (Chapter EJ) obey the char- 
acteristic power laws and show a crossover to the linear-response regime, both 
for repulsive and attractive interaction. The measured exponents agree with the 
exact values to linear order in the interaction; however the deviations are small 
up to U ~ 1.5 because we have incorporated the flow of the interaction vertex 
|aemmss 20041 Andergassen 2005| . A double barrier features several distinct power 
laws in the conductance as a function of temperature |meams 20051 Iemabms 2005j . 
We reproduce the exponents for the limiting cases treated in earlier works but in 
addition observe non-universal behavior in the intermediate parameter and temper- 
ature regions. This serves to clarify a contentious issue in the literature where two 
different universal scaling laws have been claimed for a certain parameter range: our 
data agree qualitatively with the numerically exact QMC |HuglefeEgger 2004| but 
suggest an interpretation in terms of a non-universal crossover regime. In conclu- 
sion, it is remarkable that our simple approximation captures on an equal footing 
effects which originate from very different physical processes. 

As an outlook there are several promising extensions of the current scheme, (i) The 
electron spin should be included to make the model more realistic |Andergassen 2005] . 
Preliminary results indicate that for the typical size of carbon nanotubes (N ~ 10 4 ), 
the results are even farther from the asymptotic behavior than in the spinless case, 
so the accurate treatment of the non- universal behavior is even more important. 
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(ii) X. Barnabe-Theriault studied junctions and rings pierced by a flux, and con- 
nected to several interacting leads |bsms 2UU5a| Ibsms 2005bj . He showed that the 
interpretation of the complicated behavior of such systems is greatly simplified if 
one considers not the (real) conductance but the (complex) Green function at the 
interfaces to the leads as the relevant observable: if plotted parametrically for dif- 
ferent impurity strengths and energy scales, one arrives at a simple flow diagram in 
the complex plane, (in) The investigation of non-equilibrium phenomena using a 
Keldysh variant of the fRG has been started [Jakobs 2004j . (iv) There are interesting 
Luttinger-liquid properties which appear only at two-loop order in the interaction, 
related to inelastic scattering. To take them into account one would need to include 
the frequency dependence of the interaction vertex or two-loop diagrams into the 
flow equation. 

The fRG provides a powerful tool to compute properties of Id lattice models where 
all microscopic parameters can be flexibly modeled. It captures the effects on many 
energy scales, yielding universal scaling as well as non-universal behavior. 
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A Heat equation 



I hope to develop a better intuition of the formal procedures and transformations 
involved in the functional formalism by repeating them on a simple, well-known 
example: the heat equation in one dimension. Given a temperature distribution 
uq(x) at initial time t = 0, the heat flow governed by the equation 



determines the temperature distribution u t (x) at any later time t > 0. 

In the left column below, we start with the solution of (|A.1J) in integral form, 
a convolution with the Green function, and transform it in several steps into the 
differential form. At the same time in the right column, we start with the functional 
integral definition ()2.13|) of the effective interaction V because the generating func- 
tional e~ v is the one formally most closely related to the temperature distribution 
Ut(x) at fixed time t — 1: 



d t u t (x) = d 2 x 



u t (x), 



u t=0 (x) = u (x), 



(A.l) 




x e 



measure: djj, t (x) = M ■ e « dx, J djJL t {x) = 1 





shift of variables: 




s=0 




—Vo[8$,6 v ] 



X 




V?=0 
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completing the square: M J e « 

= u (d s ) e ts2 e 
= u (d. 



y — \-ys 



aly 



As* 



e -V Q [5(p,5<p] e i<fi,C<p) e (ip,x)-(x,<f) 



s=0 







td 2 i 

e x e 



s=0 



e tA u (x) 



e Ac e ~ Vo[x 'X ] 



The solution of the heat equation is most conveniently formulated in Fourier space 
where the Laplacian A = — k 2 is diagonal. However, on the field-theory side per- 
turbation theory in the interaction Vq leads to an expansion of e~ Vo in powers of \i 
which corresponds to an expansion of Uq(x) in powers of position x, not momentum 
k. Therefore, we will try to solve the heat equation in real space to make the analogy 
clearer. On the way we shall obtain a generalization of Hermite polynomials as basis 
functions and see how they are obtained from a generating function, the Id analog 
of generating functionals. In the left column, we will first expand uq(x) in powers of 
x and then apply e tA to it; in the right column, we proceed in the opposite order, 
first applying e tA and then expanding in powers of x: 



u t {x) 



e tA u (x) 
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Hl(x). 



J o 



(A.2) 
(A.3) 

(A.4) 
(A.5) 



In the last line, we have used "rescaled Hermite polynomials" Hl(x) defined by 



Hl(x) 



e td *x\ 



- s2+xs = E sk 



(A.6) 



which are related to (usual) Hermite polynomials Hf.(x) 



H k (x) 



(2x) 



s 2 +2xs 



E^ 



k [x) s 



by 



Hl(x) 



(A.7) 



(-t) k / 2 H k (x/V^4t). 

Thus, once we have expanded the initial condition into modes Uq , (0), they evolve 
in time independently of each other according to the t dependence of H l k (x). As 
there are infinitely many modes in Uq{x) ~ e _v °' x '^, however, it turns out that other 
parametrizations are more efficient in our applications. 
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B.l Propagator 



For the frequency cutoff at T = and A < A < oo the flow equation for the self 
energy (J4.24J1 can be written as 



(B.i; 



2,2' 



In order to compute its right-hand side, one needs to invert the tridiagonal matrix 
T = G \iA) -S A , (B.2) 



where T is complex symmetric (not hermitean) with diagonal elements a, := iA + 
fi — Hifi, i = 1, . . . , N, and first off-diagonal elements bi := t — i = 1, . . . , N — 1. 

Note that Im(eij) = A > such that T is non-singular and its inverse well-defined. 

The inverse G A (iA) = T _1 is not tridiagonal but a full matrix which can be 
computed by standard methods in 0(N 2 ) time. However, for an interaction that 
does not extend beyond nearest neighbors on the lattice, only the tridiagonal part of 
G is required, which can be computed in 0(N) time, such that much larger lattices 
can be treated. We shall first explain how this is done and then present the resulting 
algorithm that can directly be incorporated into a computer program. 

Under certain assumptions (see below), a matrix can be uniquely factorized into 
a lower unit triangular matrix L, a diagonal matrix D, and an upper unit triangular 
matrix U ("LDU factorization"): T = LDU |ptvf IHHEj - For a tridiagonal matrix 
T the unit triangular matrices L and U are in fact unit bidiagonal: their matrix 
elements are unity on the diagonal, and only the first off-diagonal is nonzero. Since 
our T is symmetric we have L = U T . Thus we obtain a factorization of the form 



T = U +T D + W 



( 1 

Ut 1 



\ 



1 



V 



V 



\ (1 Ut 

1 ut 



d; 



■J 



V 



where the label "+" distinguishes this factorization from another one used below. 
The prescription to compute the elements Df and XJf is well known and can be 
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found for example in |ptvf 198fij . Starting in the upper left corner one proceeds to 
increasing row and column numbers until one arrives at the lower right corner of T: 

Dt:=a u Uf:=bi/Df, Df +1 := a i+1 - hU+ (i = 1, . . . , N - 1) . (B.3) 

This works well since in our case lm(Df) > A > 0, such that one never divides by 
zero. 

To compute the inverse G = T _1 , one could directly calculate (U + )~ 1 (D + )~ 1 (U +T )~ 1 . 
It is however easier and more accurate to find the inverse by solving the linear system 
of equations TG = 1, where 1 is the identity matrix, by "back substitution". To be 
specific, consider the i th column vector G.^ of G: 

6i = TG., = U +T (D + U + G.,) = U +T g u U + G., = (D + y l 9l (B.4) 

where ej is the i th unit vector. The first step is to solve the linear system U +T g.i = e\ 
for gj, and the second step to solve U + G.^ = (D + )~ 1 g i for To solve a tridiagonal 
linear system for one vector takes O(N) time, so solving for the full inverse matrix 
G takes 0(N 2 ) time. 

Now we shall derive an algorithm to compute the elements of gi and G.^. Begin 
with the last column i — N: U +T gN = &n can be solved from the first to the last 
row and gives g^ = &n- Next U + G.^ = (-D + ) _1 eAr can be solved starting from the 
last row, (5jv,jv = V-^iv- F rom there one can work upwards by back substitution, 
Gj t N = —Uj~Gj + i t N (j = 1, ■ ■ ■ , N — 1). For the other columns i < N, one cannot 
take the shortcut and has to solve both linear systems for gi and G.^. But it is 
now important to realize that for any column vector G.j+i, if we somehow know the 
diagonal element Gj+i^+i, the next element above the diagonal is 

G iti+1 = -U+G i+1>i+1 (i = l,...,N-l). (B.5) 

Thus, we have a prescription how to go up one row in G. Together with the symmetry 
of G, i.e., Gi t i+\ = (jt+iji, which follows from the symmetry of T, we get the first 
off-diagonal element one column to the left without solving the two linear systems in 
(|B.4[) . Hence, it is possible to compute directly the tridiagonal part of the inverse. 
However, there is another algorithm which is much more accurate for near-singular 
matrices at the end of the RG flow: the double factorization jMeurant 1992] . It does 
not rely on the symmetry of G but uses the complementary "UDL" factorization 

T = U~D~L~ = U~ D~U~ T , (B.6) 

where the matrix elements are obtained as 

D N :=a N , Ur :=bi/Df +1 , D; := a, - b,Ur (i = N - 1, . . . , 1) . (B.7) 
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We proceed as for the ldu factorization above and get 

G 1A = 1/L>r (B.8) 
= -Uf G M . (B.9) 

We can combine equations ()B.5jl and ()B.9)) to relate consecutive diagonal elements: 

G l+W = -G i>i+1 /Ut = G it i Uf/U+ = G iti Df / Df +1 . (B.10) 

Thus, we start with (|B.8J) and use the U~D~L~ decomposition to go one matrix 
element to the right in the inverse matrix, from the diagonal to the first off-diagonal 
(JB.9)) . while the L + D + U + decomposition allows to go down by one, back to the next 
diagonal element (jB.5|) . There is no need to compute the full inverse matrix. 

One can implement the algorithm without knowing the derivation by using equa- 
tions (jB.3|) and (|B.7|) - (jB.10)l . One can further eliminate the C/'s using equations 
(jB.3|) and (jB.7|) and implement the algorithm such that only the input vectors a^, 
bi and the output vectors G^i, G^j+i enter the temporary storage. This double 
factorization is numerically accurate to more than 10 significant digits (using dou- 
ble precision) even for large lattices (10 6 sites) and almost singular matrices with 
|aj| ~ 10~ 15 which appear at the end of the flow for half filling. 



B.2 Bubble 

Even the RHS of the flow equations 1)4.35)1 for the density-response vertex R A and the 
flow of the self energy for T > in equation 1)4.30)1 can be computed in O(N) time 
jEnss 2005) . In both cases the RHS of the flow has the form of a trace of a bubble 
with two vertices and two propagators, 

tr(UG 1 TG 2 ) = tv(UM) (B.ll) 

where we have defined the product 

M = G 1 TG 2 (B.12) 

with U, T tridiagonal matrices ("vertices") and Gi, G 2 inverse tridiagonal matrices 
("propagators"). In order to compute the trace in the end and also as a useful 
intermediate result in the flow equation, we need the tridiagonal part of M, i.e., Ma, 

M i)i+ i and M i+1)i . 

According to jMeurant 1992) the inverse of a tridiagonal matrix generally has the 
following structure: the upper triangle is spanned by two vectors Xi, yi, while the 
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lower triangle is spanned by two vectors Ui, vf. 

UuVik i > k I u 2k v 2 j k > j (lower triangle) 

Gi-,ik = \ xikUu i < k G 2 -kj = < x 2 jH2k k < j (upper triangle) (B.13) 
An i = k I A 2 fc k — j (diagonal) 

On the diagonal, the upper and lower triangles agree, and we define an abbreviation 
Aj := x-iUi = u(Vi for it. The four vectors have the following interpretation: Ui (mul- 
tiplied by V\) is the first column, Vi (multiplied by u^) the last row; Xi (multiplied 
by yi) the first row, and yi (multiplied by xn) the last column. Each of these vectors 
can be computed via the L + D + U + and U~ D~L~ decompositions which relate one 
row or column to the next or previous: 

— = — = -Lf (B.14) 

Ui v i+1 

^ = -Ur = -Ut (B.15) 

Consider first the diagonal elements Ma, then the off-diagonal elements Mj ii+1 and 
Mj +1 j will be slight variations of it. The tridiagonal matrix T = (a, b, c) has diagonal 
elements a^, above the diagonal hi and below the diagonal q: 



Ma — [G\-ikakG2-ki + G\-ikbkG2-k+i,i + G\-i :k +iCkG2-ki\ 

uiiV lk a k y2kX2i - A li aiA 2 i + ^ ynx lk a k u 2k v 2i \ 

k<i k>i 

+ uiiV lk b k y2,k+iX2i + ^2 yu x ikhu2,k+i v, 2i } 

k<i k>i 

+ UuVl,k+lCkV2kX2i + y^|/nXi i fc + lCfcM 2 fc^2»| 

k<i k>i 

= Qf - AiiOiAa + Qf 

where Qf are the terms on the left side inside each curly bracket (lower k indices) 
and Qf those on the right side (higher k indices), so up to now we have just split 
the terms in the k sum into two groups. The point of defining the Q's like this is 
the following: Qf +1 can be computed from Qf in 0(1) time and likewise Qf from 
Qf + i, so going from Ma to Mj+i^+i is an 0(1) operation. Thus we can compute the 
whole diagonal of M in O(N) time! 

The recursion relation for Qf is read off from the above partition: 

Qf = u u x 2 i ^2 ( v ikakV2k + v lk b k y 2>k +i + v i,k+i c kU2k) + Ai i a i A 2i 

k<i 
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Qt +1 = u lti+1 x 2 , i+ i (v lk a k y 2k + v lk b k y 2M1 + v 1M1 c k y 2k ) + A lti+ ia i+1 A 2)i+ i 

k<i+l 
Ul i+i X 2 i + i + 

Uli x 2i 

+ Ai i j + iaj + iA2,j+i 

= LiiQi U 2i — L li AubiA 2t i + i — Ai : i + iCiA 2 iU 2i + Ai^+iOj+iA^j+i 

and likewise for 

Qi = y u v 2i 22 { x ikO-kU 2 k + x lk b k u 2 , k+1 + x ljk+ ic k u 2k ) 

k>i 

Qi = UtiQ~i+\Lti - A H 6iA 2 ,i + iL 2 f i - U^A hi+1 CiA 2i + A H a;A 2 ;. 
In the same way the off-diagonal elements M iji+ i and M i+ i ti are determined: 
Mi — ^ [Gi;i k a k G 2 . k ^ i+ i + Gi-i k b k G 2 . k+ i^ + i + Gi. ijk+ ic k G 2 . kti+ i[ 

k 

= UuVi k a k y 2k x 2;i+1 + ^ yux lk a k u 2k v 2ti+ i j 

+ 1^^1^1/0^2/2^+1^2,2+1 + UiiViibiy 2 , i+1 x 2 , i+ i + ^ yiiXi k b k u 2tk+1 v 2)i+1 } 

+ U li Vl,k+lCky2kX2,i+l + Wli^l,i+lCi2/2i^2,i+l + ^ yiiXl, k+ lC k U 2k V 2 ,i+i } 

= + (A li 6 i A 2)i+1 ) + ([/+A 1)i+lCi A 2i f/ 2 7) - C/+Q- a 

and 



Mt+i,i — [Gi^+i^afcG^fcj + Gi ; j + i ) fc6fcG'2;fc+i,j + Gi^+i^+iCfcG^fci] 

= u lti+ iv lk a k y 2k x 2i + ^ 2/i,i+iXu.a fe u 2 fcf2i| 

fc<i fc>i 

+ Ui,i+iVikhy2,k+ix 2 j + u lti+1 v u biy 2t i + ix 2i + y^ j yi, i+ ix lk b k u 2M1 v 2i ^ 

k<i k>i 

+ | ^ u 1)i+ iv 1)k+ ic k y 2k x 2i + Ui.i+ifi.i+iCj^x^ + yi,i+ixi tk+ ic k u 2k v 2i j 

= -L u Qf + (Li i A li 6 i A 2ji+ iLj) + (A lii+lCi A 2i ) - QT +1 L+. 

After the full derivation, let us summarize the algorithm. We observe that certain 
combinations of terms appear frequently so we define the following abbreviations: 

Ai := A u aiA 2i (B.16) 



107 



B Efficient computation of tridiagonal loops in 0(N) 



B, t := A u kA w (B.17) 

C % := A M+lQ A 2i . (B.18) 

At the beginning of the algorithm, Qf and Qf need to be computed via the recursion 
formulas 

Qt +1 := LuQfVZ - LfB t - QU^ + A i+1 (B.19) 

Qf := U+Q^L* - B.L+ - Vtfii + M (B.20) 

with initial conditions 

Qt = A 1 (B.21) 

Q N = A N . (B.22) 

Then the tridiagonal components of M can be computed in any order: 

M M := Q+ -A + Qf = Qf - fl,Lj - U+Q + UtQf +1 L+ (B.23) 

M M+1 := B t - Q+Ui - U+Qf +1 + U+QU^ (B.24) 

M i+hi := d - L~ u Qt - Qf +1 L+ + Lf^L+i. (B.25) 

Finally, the trace (|B.11|I is 

tx(UM) := ^(%M 44 + Ui 
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